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Description 

Field of the Invention 



5 [0001] The invention generally relates to techniques for analyzing biological samples such as DNA, RNA, or protein 
samples and in particular to techniques for analyzing the output patterns of hybridized biochip microarrays. 

Background of the Invention 

10 [0002] A variety of techniques have been developed to analyze DNA or other biological samples to identify diseases, 
mutations, or other condrtions present within a patient providing the sample. Such techniques may detemiine, for 
example whether the patient has any particular disease such as cancer or AIDS, or has a predisposition toward the 
disease, or other medical conditions present in the patientDNA-based analysis may be used either as an in-vitro or 
as an in-vivo control mechanism to monitor progression of disease, assess effectiveness of therapy or be used to 

15 design dosage formulations. DNA-based analysis is used verity the presence or absence of expressed genes and 
polymorphisms. 

[0003] One particularly promising technique for analyzing biological samples uses a DNA-based microan-ay (or mi- 
croelectronics biochip) which generates a hybridization pattern representative of the characteristics of the DNA within 
the sample. Briefly, a DNA microarray includes a rectangular an-ay of immobilized single stranded DNA fragments. 
20 Each element within the an-ay includes few tens to millions of copies of identical single stranded strips of DNA containing 
specific sequences of nucleotide bases. Identical or different fragments ofDNA may be provided at each different ele- 
ment of the array. In other words, location (1.1) contains a different single stranded fragment of DNA than location (1 .2) 
which also differs from location (1 .3) etc. Certain biochip designs may replicate the nucleotide sequence in multiple cells. 
[0004] DNA-based microarrays deploy chemiluminiscence. fluorescence or electrical phenomenology to achieve the 
25 analysis. In methods that exploit fluorescence imaging, a target DNA sample to be analyzed is first separated into 
individual single stranded sequences and fragmented. Each sequence being tagged with a fluorescent mari<er mole- 
cule The nagments are applied to the microarray where each fragment binds only with complementary DNA fragments 
already embedded on the microan-ay. Fragments which do not match any of the elements of the microarray simply do 
not bind at any of the sites of the microanray and are discarded during subsequent fluidic reactions. Thus, only those 
30 microan-ay locations containing fragments thai bind complementary sequences within the target DNA sample will re- 
ceive the fluorescent molecules. Typically, a fluorescent light source is then applied to the microarray to generate a 
fluorescent image identifying which elements of the microan-ay bind to the patient DNA sample and which do not. The 
image is then analyzed to determine which specific DNA fragments were contained within the original sample and to 
detennine therefrom whether particular diseases, mutations or other conditions are present in the patient sample. 
35 [0005] For example, a particular element of the microanray may be exposed to fragments of DNA representative of 
a particular type of cancer. If that element of the array fluoresces under fluorescent illumination, then the DNA ofthe 
sample contains the DNA sequence representative of that particular type of cancer. Hence, a conclusion can be drawn 
that the patient providing the sample either already has that particular type of cancer or is pertiaps predisposed towards 
that cancer. As can be appreciated, by providing a wide variety of known DNA fragments on the microanray, the resulting 
40 fluorescent image can be analyzed to identify a wide range of conditions. 

[0006] Unfortunately, under conventional techniques, the step of analyzing the fluorescent pattern to detemriine the 
nature of any conditions characterized by the DNA is expensive, time consuming, and somewhat unreliable for all but 
a few particular conditions or diseases. One major problem with many conventional techniques is that the techniques 
have poor repeatability. Hence, if the same sample is analyzed twice using two different chips, different results are 
45 often obtained. Also, the results may vary from lab to lab. Consistent results are achieved only when the target sample 
has high concentrations of oligonucleotides of interest. Also, skilled technicians are required to prepare DNA samples, 
implement the hybridization protocol, and analyze the DNA microan-ay output possibly resulting in high costs. One 
reason that repeatability is poor is that the signatures within the digitized hybridization pattern (also known as a "dot 
spectrogram-) that are representative of mutations of interest are typically very weak and are immersed in considerable 
50 noise. Conventional techniques are not particulariy effective in extracting mutation signatures from dot spectrograms 
in low signal to noise circumstances. Circumstances wherein the signal to noise ratio is 0 to strongly negative (-2 to 
-30dB) are particulariy intractable. 

[0007] Accordingly, it would highly desirable to provide an improved method and apparatus for analyzing the output 
of the DNA microan-ay to more expediently, reliably, and inexpensively determine the presence of any medical condi- 
55 tions or concerns within the patient providing the DNA sample. It is particulariy desirable to provide a technique that 
can identify mutation signatures within dot spectrograms even in circumstance wherein the signal to noise ration is 
extremely low. It is to these ends that aspects of the invention are generally drawn. 

[0008] Refen-ing now to FIG. 1 . conventional techniques for designing DNA microarray chips and for analyzing the 
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output thereof will now be described in greater detail. Initially, at step 100. fluorescently labeled primers are prepared 
for flanking loci of genes of interest within the DNA sample. The primers are applied to the DNA sample such that the 
fluorescently labeled primers flank genes of interest. At step 1 02. the DNA sample is fragmented at the locattons where 
the fluorescently labeled primers are attached to the genes of interest to thereby produce a set ofDNA fragments, also 

5 called "oligonucleotides'* for applying to the DNA microanBy. 

[0009] In general, there are two types of DNA microarrays: passive hybridization microan-ays and active hybndization 
microarrays Under passive hybridization, oligonucleotides characterizing the DNA sample are simply applied to the 
DNA microarray where they passively attach to complementary DNA fragments embedded on the an-ay. With active 
hybridization, the DNA array is configured to extemally enhance the interaction between the fragments of the DNA 

10 samples and'the fragments embedded on the microarray using, for example, electronic techniques. Within FIG. 1 , both 
passive hybridization and active hybridization steps are illustrated in parallel. It should be understood that, cun-ently 
for any particular microan-ay, either the passive hybridization or the active hybridization steps, but not both, are typically 
employed. Refening first to passive hybridization, at step 104 a DNA microarray chip is prefabricated with oligonucle- 
otides of interest embedded or othenA^ise attached to particular elements within the microarray. At step 106. the oligo- 

t5 nucleotides of the DNA sample generated at step 1 02 are applied to the microan^y. Oligonucleotides within the sample 
that match any of the oligonucleotides embedded on the microarray passively bind with the oligonucleotides of the 
array while retaining their fluorescently labeled primers such that only those locations in the microarray having corre* 
sponding oligonucleotides within the sample receive the primers. It should be noted that each individual nucleotide 
base within the oligonucleotide sequence (with lengths ranging from 5 to 25 base pairs) can bond with up to four 

20 different nucleotides within the microan-ay. but only one oligonucleotide represents an exact match. When illuminated 
with fluorescent light, the exact matches fluoresces most effectively and the non-exact matches fluoresce considerably 
less or not at all. 

[0010] At step 108, the DNA microan-ay with the sample loaded thereon is placed within a fluidics station provided 
with chemicals to facilitate the hybridization reaction, i.e.. the chemicals facilitate the bonding of the oligonucleotide 
25 sample with con-esponding oligonucleotides within the microan-ay. At step 110, the microarray is illuminated under 
fluorescent light, pertiaps generated using an ion-argon laser, and the resulting fluorescent pattern is digitized and 
recorded. Alternately, a photograph ofthe fluorescent pattern may be taken, developed, then scanned into a computer 
to provide a digital representation of the fluorescent pattem. In any case, at step 11 2. the digitized pattem is processed 
using dedicated software programs which operate to focus the digital pattem and to subsequently quantize the pattem 

30 to yield a fluorescent intensity value for each an^y within the microarray pattem. At step 114, the resulting focused 
an^y pattem is processed using additional software programs which compute an average intensity value at each array 
location and provides for necessary normalization, color compensation and scaling. Hence, following step 114. a dig- 
itized fluorescent pattem has been produced identifying locations within the microan-ay wherein oligonucleotides from 
the DNA sample have bonded. This fluorescent pattern is refenred to herein as a "dot spectrogram", 

35 [0011] In existing biochips that actively initiate, facilitate or selectively block hybridization, a DNA microan-ay is pre- 
fabricated for active hybridization at step 1 16. At step 118. the DNA sample is applied to the active an-ay and electronic 
signals are transmitted into the anray to help facilitate bonding between the oligonucleotides of the sample and the 
oligonucleotides of the an-ay. The microanray is then placed within a fluidics station which further facilitates the bonding. 
Thereafter, at step 122, an electronic or fluorescent readout is generated from the microan^y. When electrical output 

40 signals from the biochip an-ay are used to quantify and classify the post-hybridization output, the output signal is in- 
dicative ofthe number oligonucleotide fragments bonded to each site within the array. At step 124 the electronic output 
is processed to generate a dot spectrogram similar or identical to the dot spectrogram generated using the optical 
readout technique of steps 110-1 14. Hence, regardless of whether steps 104-1 14 are perfonned or steps 116-1 24 are 
perfomied the result is a dot spectrogram representative of oligonucleotides present within the DNA sample. Here it 

45 should be noted that some conventional passive hybridization DNA microan^ys provide electronic output and some 
active hybridization microelectronic anrays provide optical readout. Thus, for at least some techniques, the output of 
step 108 is processed in accordance with steps 122 and 124. For other techniques, the output of step 120 is processed 
in accordance with steps 110-114. Again, the final results are substantially the same. i.e.. a dot spectrogram. 
[0012] At step 126. the dot spectrogram is analyzed using clustering software to generate a gene an-ay amplitude 

50 readout pattem representative of mutations of interest within the target DNA sample. In .essence, step 126 operates 
to con-elate oligonucleotides represented by the dot spectrogram with corresponding DNA mutations. Next, at step 
1 28, the resulting digital representation of the mutations of interest are processed using mapping software to detemiine 
whether the mutations are representative of particular diagnostic conditions, such as certain diseases or conditions. 
Hence, step 128 operates to perform a mutation>to-diagnostic analyses. Finally, at step 130 the diagnostic conditions 

55 detected using step 1 28 are evaluated to further determine whether or not the diagnostic, if any. can property be based 
upon the DNA sample. Classical methods such as probabilistic estimator such as minimum a posteriori (MAP) estimator, 
maximum likelihood estimator (ULE) or inferencing mechanism may be used to render a diagnostic assessment. 
[0013] As noted above, it would be desirable to provide improved techniques for analyzing the outputs for DNA 
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microarrays to more quickly, reliably and Inexpensively yield a valid diagnostic assessment. To this end. the invention 
is directed primarily to providing a sequence of steps for replacing steps 114-130 of FIG. 1. 

Summary of the Invention 

5 

[0014] In accordance with the invention, a method according to claim 1 is provided for analyzing an output pattem 
of an arrayed microdevice to identify mutations, if any. present in a biological sample applied to the blochip. In accord- 
ance with the method, a resonance pattem is generated which is representative of resonances between a stimulus 
pattem associated with a set ofknown mutations and the output pattem of the arrayed microdevice. The resonance 
10 pattem Is interpreted to yield a set of confimied mutations by comparing resonances found therein vwth predetenmined 
resonances expected for the selected set of mutations. 

[001 51 In an exemplary embodiment, the biological sample is a DNA sample and the output pattem being analyzed 
is a quantized dot spectrogram generated by a hybridized oligonucleotide microarray. The resonance pattem is gen- 
erated by iteratively processing the dot spectrogram by performing a convergent reverberation to yield a resonance 

15 pattem representative of resonances between a predetemilned set of selected Quantum Expressor Functions and the 
dot spectrogram until a predetennined degree of convergence is achieved between the resonances found in the res- 
onance pattem and resonances expected for the set of mutations. The resonance pattem is analyzed to yield a set of 
confirmed mutations by mapping the confirmed mutations to known diseases or diagnostic conditions of interest, as- 
sociated with the pre-selected set of known mutations to identify diseases, if any, indicated by the DNA sample. A 

20 diagnostic confinnation is then made by taking the identified diseases and solving in reverse for the associated Quantum 
Expressor Functions and then comparing those Quantum Expressor Functions with ones expected for the mutations 
associated with the identified disease to verify correspondence. If no con-espondence is found, a new sub-set of known 
mutations are selected and the steps are repeated to detemiine whether any of the new set of mutations are present 
in the sample. 

25 [001 6] In the exemplary embodiment the set of nonlinear Quantum Expressor Functions are generated are follows. 
A set of mutation signatures representative of the pre-selected set of known mutations is input. A representation of a 
microan-ay oligonucleotide pattem layout for the microarray. from which the dot spectrogram is generated, is also input. 
Then a set of resonant interaction parameters are generated which are representative of mutation pattem interactions 
between elements of the mlcroan^y including interactions from a group including element-to-element interactions. 

30 element-to-ensemble interactions, ensemble-to-element interactions, and ensemble-to-ensemble interactions. Then 
the set of nonlinear Quantum Expressor Functions are generated from the set of resonant interaction pattems by 
matching selected hamionics of the power spectral density (PSD) amplitude of a coded mutation signature, conre- 
sponding to the pre-selected mutation set of interest, to that of a pre-determined quantum-mechanical Hamiltonian 
system so that stochastic and detenninistic time scales match, and the time scales couple back to noise statistics and 

35 degree of asymmetry. 

[001 7] Also in the exemplary embodiment, the dot spectrogram is differentially enhanced prior to the generation ofthe 
resonance pattem by refocusing the dot spectrogram to yield a re-focused dot spectrogram; cross-con-elating the ref- 
focused dot spectrogram; applying a local maxima filter to the conrelated re-focused dot spectrogram to yield a maxi- 
mized dot spectrogram; re-scaling the maximized dot spectrogram to yield a unifonmly re-scaled dot spectrogram; and 
40 then re-scaling the uniformly re-scaled dot spectrogram to amplifying local boundaries therein to yield a globally re- 
scaled dot spectrogram. 

[0018] By exploiting a resonant interaction, mutation signatures may be identified within a dot spectrogram even in 
circumstances involving low signal to noise ratios or, in some cases, negative signal to noise ratios. By permitting the 
mutation signatures to be identified in such circumstances, the reliability of dot spectrogram analysis is thereby greatly 

45 enhanced. With an increase in reliability, costs associated with perfomning the analysis are decreased, in part, because 
there is less of a requirement for skilled technicians. Other advantages of the invention arise as well. 
[0019] Among other applications, principles of the invention are applicable to the analysis of various arrayed biomo- 
lecular. ionic, bioelectronic. biochemical, optoelectronic, radio frequency (RF) and electronic microdevices. Principles 
ofthe invention are particulariy applicable to mutation expression analysis at ultra-low concentrations using ultra-high 

50 density passive and/or active hybridization DNA-based microarrays. Techniques implemented in accordance with the 
invention are generally independent of the physical method employed to accumulate initial amplitude information from 
the bio-chip array, such as fluorescence labeling, charge clustering, phase shift integration and tracer imaging. Also, 
principles of the invention are applicable to optical, optoelectronic, and electronic readout of hybridization amplitude 
pattems. Furthemiore. principles of the invention are applicable to molecular expression analysis at all levels of ab- 

55 straction: namely DNA expression analysis, RNA expression analysis, protein interactions and protein - DNA interac- 
tions for medical diagnosis at the molecular level. 
[0020] Apparatus embodiment according to claim 20 is also provided. 



4 



EP 1 145 181 B1 

Brief Description of the Drawings 

[0021] The features, objects, and advantages of the present invention will become more apparent from the detailed 
description set forth below when taken in conjunction vwth the drawings in which like reference characters identify 
5 correspondingly throughout and wherein: 

Fig. 1 is a flow chart illustrating conventional passh^e and active hybridization DN A microan-ay analysis techniques. 
Fig. 2 is a flow chart illustrating an exemplary method for analyzing the output of a hybridized DNA microarray in 
accordance with the invention. 
10 Fig. 3 graphically illustrates the method of Fig. 2. 

Fig. 4 is a flow chart illustrating an exemplary method for generating Quantum Expressor Functions for use with 
the method of Fig. 2. 

Fig. 5 is a flow chart illustrating an exemplary method for preconditioning the output of a hybridized DNA microarray 
for use with the method of Fig. 2. 
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Detailed Description of Exemplary Embodiments 



[0022] With reference to the remaining figures, exemplary embodiments ofthe invention will now be described. The 
invention will be described primarily with respect to an exemplary method for analyzing mutations signatures within 
20 output pattems of hybridized microarrays generated using DNA samples, but principles of the invention apply to the 
analysis of other protein-based samples or to other types of output pattems as well. 

Overview 

25 [0023] Briefly, the exemplary method exploits, among other features: (a) a novel representation. Interpretation and 
mathematical model for the immobilized oligonucleotide hybridization patterns, represented via a dot spectrogram; (b) 
a new "active" biomolecular target detection and discrimination method based on quantum resonance interferometry. 
and (c) a new spatial hashing function that yields accurate diagnostic assessment. 

[0024] To this end, the exemplary method exploits a fundamentally different computational paradigm for mutation 
30 expression detection in pre-enhanced dot spectrogram realizations. The method is an innovative modification to dy- 
namically anayed quantum stochastic resonance (QSR) for discrete system analysis. The an-aying strategy is a function 
of the expression pathway of interest. The method depends on the molecular diagnostic spectrum being addressed. 
Coupled quantum resonators are employed to signiflcantly enhance signal-to-noise (SNR) performance and fuse mul- 
tiple synthetic renormalized dot spectrogram realizations to better detect prespecified biomolecular expression pat- 
35 terns. 

[0025] Moreover the exemplary method exploits an enhancement in previous extensrons to classical stochastic res- 
onance (SR) and array enhanced SR (AESR) in signal processing and sensor data analysis. Stochastic resonance is 
a phenomenon wherein the response to a sensor, modeled in terms of a bistable nonlinear dynamical system, is en- 
hanced by applying a random noise element and a periodic sinusoidal forcing function. SR occurs when the SNR 

40 passes through a maximum as the noise level is increased. 

[0026] Thus an important aspect of the exemplary method involves the coupling of transformed and preconditioned 
discrete microanay outputs to a mathematical model for a quantum-mechanical dynamical system with specific prop- 
erties. When driven in a particular manner, the coupled system exhibits a nonlinear response that conresponds to 
detection of phenomena of interest. The method exploits modulation of observables from a "base" (canonical contin- 

45 uous dynamical system), so that a selected set of spectral properties match a similar selected spectral properties of 
a discrete spatial tessellation substructure from an amplitude spectrogram derived from bioelectronic observables. The 
method further exploits the concept of convolving a discrete spatial system (derived from base mutants of interest) 
with a continuous asymmetric temporal system to derive a spatiotemporal input to further convolve with another discrete 
spatial projection (of an inherently partially stabilized spatiotemporal system). 

50 [0027] Hence key components of the exemplary biomolecular detection method are: (i) selection of a basis system; 
(ii) generation of designer Quantum Expressor Function (QEF) for coupling with the substrate to be analyzed; (iii) 
generation of a Hamiltonian to describe relaxation dynamics of the coupled system; (iv) modulation of resonance pa- 
rameters to enforce eariy resonance; (v) and exploitation of resonance suppressors to verify detection. 

55 The Method 

[0028] Refening first to FIG. 2. initially at step 200. a set of mutations of interest is selected. The mutations, for 
example, may be mutations relevant to cancer. AIDS, or other diseases or medical conditions. At step 201. precondh 
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tioner transforms are generated based upon the selected set of mutations. The preconditioner transfbmis are provided 
to convert mutation nucleotide sequences into expected amplitude pattems in the prespecified microarray represen- 
tation, given a particular biochip layout. At step 202, Quantum Expressor Functions are generated based upon the 
Hamiltonian of a pre-selected basis system. The Quantum Expressor Functions are designed to couple the Hamiltonian 
5 for the selected basis system to a predetermined DNA microarray configuration to permit a resonance interaction 
involving the output of the DNA microarray. Resonance stimulus is generated, at step 204. using the Quantum Ex- 
pressor functions. 

[00291 What has been summarized thus far are preliminary steps performed off-line for setting up the Quantum 
Expressor Functions and the corresponding resonance stimulus. These steps need be performed only once for a given 

10 set of mutations and for a given DNA microarray configuration. Thereafter, any number of output pattems from the 
DNA microan^y may be processed using the Quantum Expressor Functions to identify whether any of the mutations 
of the pre-selected set of mutations are found therein. Preferably. Quantum Expressor Functions are pre-^enerated 
for a large set of mutations and for a large set of DNA microan-ay pattems such that, for each new DNA microan-ay 
output pattern from each new patient sample, the presence ofany of the mutations can be quickly identified using the 

15 predetermined set of Quantum Expressor Functions. In general, the aforementioned steps need be repeated only to 
update the Quantum Expressor Functions to accommodate new and different DNA microanray pattems or to if new 
mutations of interest need to be considered. 

[0030] At step 206. an output pattern (refen-ed to herein as a Dot Spectrogram) is generated using a DNA microan-ay 
for which Quantum Expressor Functions have already been generated. At step 208, the dot spectrogram is precondi- 

20 tioned to yield a dot spectrogram tesselation (DST) to pennit exploitation of a resonance interaction between the dot 
spectrogram and the Quantum Expressor Functions. The actual resonant interaction, which involves convergent re- 
verberations, is performed at step 210 until a pre-detenmined degree of convergence is achieved. Once convergence 
is achieved, a resulting resonance pattem is processed at step 212 to identify any mutations represented thereby. As 
will be described below, step 212 is rendered trivial by virtue of the aforementioned resonant interaction which is based 

25 upon Quantum Expressor Function already conrelated with the pre-selected mutations. Hence, no complicated analysis 
is required to interpret the resonance pattem to identify the mutations. Next, at step 214, the mutations are mapped 
to corresponding diseases and conditions to thereby identify any diseases or conditions that the patient providing the 
sample being analyzed is afflicted with. Again, this is a fairly trivial step. Finally, at step 216, diagnostic confimiation 
is preformed to verify that the diseases or conditions are present in the sample. This is achieved by starting with the 

30 found diseases or conditions and then performing the steps of the method in reverse. 

[0031] FIG. 3 graphically illustrates the operation ofthe method of FIG. 2 whereby a Quantum Expressor Function 
300 is generated based on a mutation set 302. A dot spectrogram 304 is applied to the Quantum Expressor Function 
via an interferometric resonance interaction 306 yielding a resonance pattern (not shown) from which mutations sig- 
natures 308 (representative of mutations present in the sample from which the dot spectrogram was generated) may 

35 be identified by comparison with mutation set 302. 

[0032] Now details of the various steps of the method will be provided. 

Mutations Sets 

40 [0033] Refening again to FIG. 2, the mutation set of interest generated at step 200 is selected by identifying oligo- 
nucleotides representative ofthe {Z) mutations of interest. Each oligonucleotide is represented by v(l,j) which is given 
by [ooai ... 0^1, where a = {A,C,T,G} base associated with each an-ay cell [a.b] where ^0<k< 25. The entire set of 
unique oligonucleotides denoting mutations of interest, A(l.m). is given by [5o5i. ... where 5 = {A,C,T.G} length | A | 
= length |\|/|. and 0<|| A - v || < k, and the designed in y^m) oligonucleotide sequence is a perfect complement to only 

45 A(l,m)foralll.m. 

[0034] As part of step 200, an oligonucleotide table is generated which contains the oligonucleotide sequences as- 
sociated with each mutation of interest identified by row and column location (ij). The oligonucleotide table is provided 
for subsequent use at step 212 to map locations within the dot spectrogram wherein resonance occurs at step 210 to 
oligonucleotides such that mutations present in a sample being analyzed are easily identified. Also as part of step 200, 
50 a mutation table is generated which contains the diseases associated with each mutation of interest. The mutation 
table is provided for subsequent use at step 214 to map mutations identified at step 212 to specific diseases or other 
medical conditions such that the diseases can be easily identified. 



Preconditioner Transforms 



55 



[0035] Given the microan-ay layout and the mutations of interest (in the coordinate system specified by the prespec- 
ified exemplary microanay layout) preconditioner transfomis are employed at step 201 to generate an amplitude pattem 
and associated point spread function at a signal-to-noise ratio (SNR)=0. 
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[0036] In the exemplary embodiment, the DNA microarray is an N by M DNA chip array wherein an element of the 

anray is referred to herein as an "oxel o(ij). 

[0037] The pre-hybridization microarray (PEBC) is expressed as: 



PEBC" ££o(U). 

1 I 



where N and M refer to the linear (row and column) dimensions of the 2-D microan^y. 
10 [0038] An inverse Dirichlet Tessellation (IDT) on the PEBC is applied such that a singular value decomposition (SVD) 
on the resulting an-ay yields the location of the mutation of interest at the SNR=0 condition. Conceptually the process 
corresponds to setting up a concave diffusion front at the mutation-centered oxel. 
[0039] The equations used are 



Xf =k^(Dtf 



where 4/7. 

A/^=/f^L(DO"'^' 

where aN= 3/7. 



[0040] The constants k^, ka and k^ are respectively set to 0.856, 0.68. and 1.34, and for biochip dimensions N, 
M>100, t is typically set to 10^ * a^/D where a = inter-oxel distance and D is a coefficient for percolation of the point 
spread! D is chosen as 0.001. These equations impose a specific IDT. i.e., rapid diffusion on the mutation-centered 
oxel, and synthetically tesseltate the biochip. 

Basis System for Quantum Expressor Functions 

[0041] To generate Quantum Expressor Functions (QEF) at step 202, a basis system for the QEF is first selected. 
The selection of the basis system and the generation of the QEF's based thereon depends, in part, and the charac- 
40 teristics of the DNA microarray. 

[0042] The numeric value associated with each oxel is given by: 

6{iJ) = a'"^ + a^.^ .4"'^ +...+av4'+ao.4° 

45 

where [a]=IA | C | T | G] take the values [0 1 1 | 2 | 3 ] respectively 
[0043] An element of the dot spectrogram is refenred to herein as a hixel: h(ij). 

[0044] The complete spatial randomness (CSR) hypothesis for the dot spectrogram spatial point pattem is violated 
by the microarray oxel layout to maximize posteriori detectability. This means that (i) Conrect Binding (CB)-intensity i. 
50 e. , number of spatial events (hixels with con-ect bindings) do not follow a Poisson Distribution; and (ii) sampled bindings 
do not represent independent random sample from a unifonm distribution on POBC. 

[0045] Compliance to CSR hypothesis would imply constant CB-intensity over the DST, as well as no interactions 
(inhibitions or reinforcement) among the non-zero intensity hixels. Both imply an idealized microan-ay design with pos- 
itive SNR perfomiance. Thus proclivity for aggregation during hybridization is assumed. Computationally, this method 
55 assumes that a 1 000 X 1 000 microan-ay (i.e. . N = 1 ,000,000 oxels) is reduced to <0(1 00) nodes with the preconditioning 
method presented below. DST algorithm on preconditioned hixel array, such as variants of Dirichlet tessellation would 
further reduce it to ~-0( 10) with no coding. Consequently, this method must be applied to ~O(10) ensembles. For 
ultra-sensitive detection, it is assumed that the ensemble basis system is degenerate, i.e.. it has (i) a bounded intensity 
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function: (ii) a bounded radial distribution function, (iii) and is anisotropic. The method is however valid if hixel intensity 
distribution can be approximated by lattice-based processes, Cox processes, Mari^ov point processes as well as ho- 
mogeneous and non-homogeneous Poisson processes. Accommodating the complete spatial randomness (CSR) con- 
ditions, within a tessellation, to analyze "degenerate, contagious distributions" is actually paradoxical to the lineage 
5 from QSR and violate temporal linearity and contiguity assumptions for amplification of certain spatial effects. 

Quantum Expressor Functions 

[0046] As noted, the Quantum Expressor Functions (QEPs) generated at step 202 are based upon the DNA chip 
10 used to generate the dot spectrogram and based upon the mutation set of interest. More specifically, as shown in FIG. 
4. the QEF is generated based upon the spin Boson basis system by first calculating the Hamiltonian for the system 
at step 402 then, at step 404, calculating harmonic amplitudes |P J for the Hamiltonian. Next, at step 406. an order 
function (OF) is generated, then at step 408, entrainment states are measured of the OF of the ground truth. Finally, 
at step 410, the OF of ground truth is modulated to yield the QEF. 
15 [00471 The use of the spin- boson Hamiltonian permits the exploitation of quantum stochastic resonance (QSR) 
phenomena. Classical QSR is an archetypal example of a phenomenon where quantum noise is exploited to drive 
order in a quantum-mechanical system, as opposed to Gaussian noise in a classical system. The appearance of a 
resonance requires an asymmetry in the energies of the two states. A rate equation can be constmcted for the system, 
such that the dynamics can be characterized in temns of transition rates 0+ and O- between the two asymmetric 
20 quantum superposition states, and when the drive frequency and the intenwell transition rates are much slower than 
the intrawell relaxation rates. The signal to noise ration (SNR) of such a superposition system is given by: 



25 SNR = 



30 



35 



55 



4 l+expl^o^'^a^o 



10048] Where Kb is the Boltzmann's constant, and T is the temperature. The sinusoidally 

e = Eq + 5e X COSCD^^ 

modulated asymmetry energy e is given by: 

C(t) = <n^ (f+ t| q^.t) n^ (t|qo,oo)> 

[0049] In the above expression the power spectmm S{(o) represents the Fourier transfomn of C(t), containing a 
40 roughly Lorentzian broadband noise background and S-function peaks at © = 0, the driving frequency CQg. and Its har- 
monics. The measured correlation function for the quantum noise is given by C(t) = <n+| (t)n+i (t+ t)>. where each n+i 
= 0 or 1, and the <, indicates an average over t over many data points i taken at equal intervals. 
[0050] The probability of being in the + quantum state at t after being in the state qo at is given by 

n^(/|q^.-oo). 

QSR occurs for an asymmetric well, but not for a symmetric energy well. 

[0051] A custom Hamiltonian. for use with step 402, which couples the above system to an ensemble of hamionic 
50 oscillators is given by 

2 ^ n 1 



[0052] The dissipative tunneling of the quantum system can be described by the above Hamiltonian. In the above 
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expression, e denotes the asymmetric energy. A is the tunneling matrix element, and o, are the Pauli spin matrices, 
and C is a harmonic oscillator creation operator with frequencies to^. 

[0053] The information aljout the effects of the environment is contained in the spectral density 

y(a»)=(>r/2)x Yn,y^S («»-<»,). 

100541 However bulk of QSR results to date are limited to either experimental obsen/ations in esoteric quantum- 

10 mechanical systems (e.g.. coherent motion of a wave packet in double quantum-well (DQW) semiconductor heter- 
ostmctures) under specific conditions or explanations of experiments under carefully controlled conditions. It does not 
appear that any algorithmic implementations or methods have been designed or implemented to enforce nonhnear 
resonance in spatial, temporal or spatio-temporal phenomena in discrete systems by coupling such systems with an- 
other mathematical system. . . ^ , .i. . . 

IS [00551 Thus, as noted above, an important aspect of the exemplary method of the invention is to couple the trans- 
formed and preconditioned discrete microarray output to a mathematical model for a quantum-mechanical dynamical 
system with specific properties. Specific exemplary parameters for use in calculating the Hamiltonian are those pro- 
r»sed by A J Legetl et al.. Reviews of Modern Physics, 59, 1 . 1 987 and A.O. Caldiera and A. J.Legett Annals of Physics 
149 374 1983 The parameters are important only for an offline simulation of this spin Boson system on a digital 

20 comW. The empirical obseivables are then collected and used to estimate and compute spectral properties, which 
are actually used by the method. . .„ , ... 

[0056] It is believed that any experimentally or analytical valid parameters for the above system will work with the 
technique because the robustness ofthe method depends only on the bulk and qualitative property shown by ttiis 
system and the properties of its power spectra. This is an important point since the method is actually substantially 

25 immune to nuances and specifies of the actual driving mathematical system. 

[0057] The haimonic amplitudes determine the weights of 5 spikes of an averaged spectral power density in an 
asymptotic state Soo(co). e refers to the coupling strength and defines the equilibrium state in the absence of dnving 

[OOOT] The power amplitudes ti„ in the mth frequency component of asymptotic state space are calculated at step 
30 404 using 

n„(n.e) = 4n [PJCiX^ 

35 and the phase shift is given by 



40 



45 



,^ , . ..r im(p,(n.r)) 1 



[0059] The analytic for the external force is given by 



50 [00601 The parameters y, eo are predetemiined and are design specific. Typically, values of 0.001 and 0.0001 are 
used for Y and eo respectively. In the above expression |PJ determine the weights of the 8 spikes of the averaged 

spectral power density. ,,^«^r%\ i-» j« 

[0061] For particular applications, the QEF is designed by matching the power spectral density (PSD) amplihjde of 
coded mutation signature to that of the spin-boson system described above so that stochastic and detenministic time 
55 scales match and so that the time scales couple back to noise statistics and degree of asymmetry. The method employs 
a fully automated iterative conjugate gradient relaxation method for Spectral matching between asymmetric base sys- 
tem and coded mutation signature. The actual detemiination of the QEF depends on the specifics of bioelectronics 
substrates used for actual analysis. The method is however generalizable to all or almost all arrayed embodiments. In 
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addition, the method is highly scalable to array dimensions (as the offline design trade-space time does "Ot matter to 
computational complexity). Since the system is an overdetemiined coupled system, convergence cntena and stability 
of relaxation method directly relates to downstream resonance effectiveness. 



Order Function 



[00621 At step 306. the order function (OF) of ground truth is generated as follows. The order function (OF) is for 
ground truth wherein ground truth represents a state wherein a positive signal to noise ratio (SNR) is expected for hixel 



intensities of selected oxels. 
[0063] A general form for the classical Order Function (OF) is 



where i = and 



and 0<e<1. 

[0064] The following equation 

0 

is referred to as the Diado integral. . ^ . 

100651 The OF is calculated using an oixler match which implies the variance of the density function tor a specific 

exponential family approaches 0 (within 0.000001 - 0.0001). The free energy for the density function is given by 



Where and represent state vectors. represents the random observable in symmetric bilinear form, and 
X denotes the characteristic function. . ^ ^iDu..^ho*«n 

[00661 The detailed form and the specific results for solving the system are given in D.C.Brody and L.R Hughstoru 
Geometry of Quantum Statistical Inference". Physics Review Letters. 77(14). pp. 2851-2855. 1996. An important aspect 
of the exemplary method of the invention in using the result to derive the QEF which gives unique advantages in the 
enhancement process. 

Entrain ment States 

[00671 Thus the OF derivation is based in Diado's theory of multibranch entrainment of coupled nonlinear oscillators 
wherein a number of different entrained states co-exist. At ground truth POBC locations, an a priori measurement of 
entrainment is perfomied (step 408). This is done by approximate Daido Integral () with Z, = PSD maxima at regjons 
where DST boundary > desired detection threshold. (For more infomiation regarding the Daido Integral, see Physics 
Review Letters. 77(7).1406-1411 .) Premultlpller constants are used to ensure that meets the following tests: 

maxima power spectrum density (PSD) matches; and 
. L2 norm on even hannonics < e where e = NS-MRF Barrier assuming GSR assumptions. 

[0068] The notion of entrainment states is exploited, in part, because the method treats the hixel dot spectrogram 
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as a special case of coupled nonlinear oscillators in equilibrium. However, due to device imperfections, hybndization 
degradation and other limitations number of entrained binding states (i.e. . inconrectly hybridized) coexist. A PSD match 
is desired where the OF is Z^peaked (but single peaked around each MRC-hixei). Absence of a single peak implies 
perfect or lack of hybridization. It also defines the resonance loci for this method (i.e.. where maximum SNR enhance- 
5 ment) is obtained. 

Ground Truth Modulation 

[00691 The OF of ground truth is modulated at step 410 to yield the QEF as follows. Under controlled calibration, as 
10 stated above, maximal SNR enhancement (optimal resonance) is achieved when OF yields a single peak. It is a im- 
portant design point for matching PSD of coupling spin Boson system to the synthetic QEF. The specific form of the 
QEF to be used is the generic OF shown above. So the exemplary method exploits two connotations of OF: (a) para- 
metric fomi for the QEF (that is closer to the classical forni) and (b) as exponential attractor for a dissipative system. 
The two OF's are then recoupled and convolved with the preconditioned dot spectrogram (see below). 
^5 [0070] The resulting QEF generated via the steps of FIG. 4 is given by: 

Typically, only first three PSD peaks are considered for spectral matching. 

[0071] Preferably the QEF is represented digitally using a matrix or an^y having the same number of elements as 
25 the dot spectrogram to be analyzed. 

Phase Space Resonance Stimulus Using Palm Operators 

[0072] Refemng again to Fig. 2, if the dot spectrogram generated at step 206 is not a phase-space representation 
30 of the output of the hybridization chip, then it is desirable to convert the preconditioned dot spectrogram (generated at 
step 208) and the QEF into phase space to facilitate a phase space resonance interaction. In other embodiments, 
though an amplitude-based resonance interaction is perfomied and hence it is not necessary to convert to phase 
space 'in still other implementations, other types of resonance interactions may be employed. In the following, it is 
assumed that the dot spectrogram has phase space components and hence the following conversion steps are applied 
35 to the QEF at step 204 and similar steps are applied to the dot spectrogram following image preconditioning. 

[0073] A phase embedding operator, r, is applied to the hixels corresponding to the coded base mutation set such 
that hixel values now con-espond to angles and not to intensities. These values are cyclic, with absolute magnitude of 
the phase image having no meaning. The relative magnitudes are more significant. So if X is a phase image and a is 
any constant then X+a (mod color_map_scale) is a valid descriptor for the phase image. More importantly, the differ- 
40 ence between a max_hixeLinlensity and min_hixetjntensity is 1 and not color_map_scale. The difference between 
phase values at two hixels i and j is X(i)-XG) (mod cotor_map_scale). . , • ^ 

[0074] The phase embedding operator is designed such that transitive closure between any two hixels is maintained, 
i.e.. there is an accumulated phase function v for which 



45 



50 



X(i) = y(j) (mod color_map_scale) 



[0075] The phase values of v accumulate rather than cycling back to 0. Note that y can be approximated at least 
locally by a linear function f with deviation error function en'(p) having values relatively dose to zero for which 



V(P) = ^(P) + err(p) 



[0076] This function f is the analog of an amplitude image having a constant intensity and en-(.) is the analog of the 
deviation from this constant amplitude. It is assumed that err(p) never crosses a phase discontinuity and can be treated 
like a real-valued function. If values are computed for this deviation function en-(.). then len-(p)I would specify an intensity 
image whose local maxima might indicate detection of expressed oligonucleotides of interest. 
[0077] Each phase angle con-esponds to a point on the unit circle in the coordinate lane. This is a one-tOK)ne mapping 
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10 



15 



20 



25 



?Srwm?Hnets token which is perpendicular to the vector v and passes through the point v. then mteWne wi« 
quadrant or a 90 degree arc: 

|v| = cos(45 degrees) = 0.7071 
,v| is referred to herein as the average magriitude arid the minimum a^^^^^^ 

S) anc* e = F-L, the average of e over should be N(p.m.n) to be close to 0. Hence, approx.mately. 



ave(F,p.m,n) = ave(L+e,p,m,n) 
30 = ave{L,p,m,n) + ave(e,p,m,n) 

= L(p) + ave(e,p,m,n) 
= L(P) 

35 

10081] Since t(p) = F(p) (mod), the average values of these two functions over N(p.m.n) should also satisfy this 
relationship and hence, approximately, 

40 ave(f,p,m,n) = ave(F,p,m,n) (mod) 

= L(p) (mod) 

45 10082] in other words, it is reasonably expected that ave(f.p.m.n) recovers an approximation to the underiying line- 
arity. If this is valid, then 

|ave(f,p.m,n)-f{p)| 



50 



55 



is a measure of the deviation of f at p from this undertying linearity and these values form an intensity image 

cycles. The remaining sections presents alternatives based on this assumption. 
[0084] Define the first vertical difference of a phase image flo be 
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dv(f.p) = f(i+j)-f(ij)(mod) 
where p=(i j). Similarty. define the first horizontal difference of a phase image f to t>e 

5 

dh{f.p) = f(i.j+1)-f(iJ)(mod) 

[00851 These both define new phase Images derived from f and are discrete analogs to the first derivative. The range 
10 of values for dv and df are also from -180 to + 180. ^ ^ c 

[00861 Let F be the accumulated phase function corresponding to f and let L+e=F be as described above. For the 
sake of convenience, assume that p is translated to the origin so that p=(0.0). The linear function L is expressed as 



15 



40 



45 



L(ij)=Vi+H*j+K 



for some constants V. H, and K. IfL is a reasonable approximation to F. then it is expected that the first vertical 
difference of f to be an approximation to the slope V and the first horizontal difference of to be an approximation to the 
slope H. Hence it is possible to estimate V and H by averaging dv and dh over some appropriate rectangle centered 
20 at p. Let 

V = ave( dv(f,p). p.m.n) 
H* = ave( dh(f,p). p.m.n) 

be these estimates. Note that these are phase averages. Probably m and n should be chosen so that m/n is 
roughly equal to H'A/'. This ratio would imply that the total phase accumulation across the rectangle N(p.m.n) is about 
30 the same in both the vertical and horizontal directions. For example, the rate of change of phase along the i-axis 
(vertical) may be twice the rate of change along the j-axis (horizontal). In this case, a 5x10 rectangle of hixels may turn 
out to be the most appropriate neighbortiood. 

[0087] If a good estimate K' is computed for K. then an approximation to L is: 

35 L'(i,j)=V'*i+H'*j+K' 

[0088] If L* is a good enough approximation to L and F is linear enough, then the accumulated phase function F is 
reconstructed as follows: 

e'(P) = f(P)-L'(P) (mod) 

and 

F(p) = L'(p) + e'(p) 



[0089] Note that this may require that the error function e' never deviate far enough from 0 to cross a phase discon- 

[0090] With an accumulated phase function F now explicitly computed on some rectangle N(p.m.n). standard nu- 
merical analysis techniques are applied, such as linear regression, to find a function L- which best approximates F 
over this rectangle. The value 

dev(f.p,m.n) = |F{p) - L-(p)| 

can be taken as a measure of the deviation at p of the image f from its "expected" value at p. In particular, this 
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defines a new image whose values are intensities rather than phase angles and whose local maxima ran be analyzed^ 

foMiT TO ~S a good estimate K for K in order to detemilne the function U choose 

a phaie value between -180 and *180 which minimizes the average value over the rectangle N(p.m,n) of either 



|f(P)-L'(P) (mod)l 



(f(P>-L'(P)('nod))* 



100921 If the accumulated phase function F corresponding to the phase image f is approximately a linear function 
and asalsi the fir^tdifferencesare approximately o^nslant. then th^ 

15 be approximately 0. 

[0093] Define the second vertical difference of a phase image f to be 

d2v(f.p) = lf(i+1 j) - f(i j) (mod)l - [f(i j) - f(i-1 j) (mod)] (mod) 

and the second horizontal difference of a phase image f to be 

d2h(f.p) = [f(ij+1)- f(i j) (mod)]- If(i.i)- f(i,j-1) (mod)] (mod) 

[00941 This corresponds to an indexing of a one dimensional difference table that looks like the following: 



20 



30 



35 



45 



50 



55 



aO al a2 a3 a4 a5 a6 
bO bl b2 b3 b4 b5 
cl c2 c3 c4 c5 

where bO=a1-aO. b1=a2-a1 .... and c1 =b1-bO, c2=b2-b1 
[0095] Let F=L*e be define as earlier. If e were identically 0. then F would be linear and both second differences d2v 
Lnd d2h would also be identically 0. If e is not identically 0. the second differences w,ll necessanly f o devjate from 
0 however, not necessarily in the same place as e. If |e| attains a local maximum at a h.xel p and h« .s assoaated 
with a very localized deviation from 0. as would be expected of a signal from a mutation of interest, then d2v and d2h 
should also deviate from 0 fairly nearby to p. ^ ^ . _ hiv<je 

[00961 The deviation in d2v and d2h associated with a deviation in le| may be spread out over a ""mter o^ix^s^ 
TOs will not only make If harder to detect a significant change in d2v and d2h but also make it ha^er to locate^ 
I^nteTof t^ s change using a local maxima finding filter. It may be possible to compensate for this by ^^9^ 
S2h or |d2v| and |d2h| over rectangular regions N(p,m.n). This should tend to sharpen the peaks d2vand d2h^ 
Correiatt^ several dikrent images generated using d2vandd2h may alsohelpwith^^^^^ 

algorithm will automatically average nearby local maxima. , ^ ^ .^rSK^n Ronoatahiiitv of the 

[00971 Thus a method for achieving induced coupling for phase space has been descnbed. ^^Pf^'^^''^"' 
ov^ l method hinges significantly on the aforementtoned coupling as It induces analytic redundancy va complex 
W^aMon a't conLnt feed fon«a,d computational cost. Hence, the coupling is advantageous for at least Nvo 
re~ 1) it applies to any microarray devfee that may provide a phase or cyclical (modulo) input as oPPosed to 
ar^Xde inpuTand 2) it is important L use with active hybridization devices which will have an elementof bu^t in 
^nC mat ^11 have a phase representation. The coupling is actually introduced into the exemplary method in three 
possible ways. 

(a) iflhe entire resonance computation proceeds in the phase space as opposed to amplitude space. Then it is 
necessary to transfomi QEF and post-hybridization microarray output to phase space. The coupling matenal pro- 
vides the transformations to do so. 
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(b) if the microarray output is already in the phase space then it will only be applied to convert the Q^R 

(c) if different oxels are designed to solve for different aspects of the analysis, the coupling permits a combination 
of outputs to reach a final conclusion. 

[0098] Entrained states are employed, in part, to precisely compensate for that pointy So ^ireadV in^^ ^nethod 
synthetic decomposition and coupling are provided. But it is used to reject spunous candidates as opposed to Hight 



up 



more oxels with potential match. 



[00991 Again the step of converting to phase space nray be optional depending upon the implementation and is 
applied, if at all. after preconditioning and/or before introducing the QEF to the resonant interadion s ep^ 
,0 [0100] Also, a major limiting restriction in QSR that is avoided by the exemplanr method pertains to "latchjng the 
Stochastic and deteirninistic time scales in "domain system" and the external coupling a^metnc dynam-cal^em 
since this has limited applicability to continuous data. By replacing the time scale match requirement with ensemble 
spa^a J.atis«cs (stn^Cure of oligonucleotide itself) expressed via the afo-Bmentioned generators an entire^n«^^^^ 
oLalysisfordiUtesystemsisenabledByderiving the condition thatlabeledc^i^^^^^ 
15 to maximize QEF generator capacity, a properly constmcted and conditioned QEF will be able to extract and enhance 
an entire dass of mutations over and beyond single mutation detectabihty. , . . . • k.. 

[0101] The resulting phase-space representation of the resonance stimulus generated at step 204 is given by. 



so 



55 



= (D2h{QEFMRc i) ®n»d « D2v(QEF^,rcj) ^ QEF subanay elements. 



l>iesonance stimulus " ^^"'"V^'^'^MRC. 



(0102] AS with the QEF. the resonance stimulus pattern is preferably represented digitally using a matrix or array 
having the same number of elements as the dot spectrogram to be analyzed. 



25 Generation of The Dot Spectrogram 



[0103] Continuing with the detailed description of the steps of FIG. 2, as noted a dot spectrogram is generated at 
step 206 for a sample from an N by M DNA chip array wherein an element of the array is an ofy/ A 6-a 

manufacturing process accuracy in microarray design is assumed. Each array cell amplitude is given by *0.l) for.. 1 
to N and i- 1 to M. Let denote the a priori known oligonucleotide given by looai-ak). where a - {A,C,T.e) t>ase 
associated with each array cell [a.b] where 10 < * < 25. The complimentary strand, derived from unknown sample is 

foSr ^[.ri^slhybridization microarray is treated mathematically using the machinery of equations with aftereffed 
Earhixel giv^ by *(i.j) is represented as a duster of dynamical systems of PO'e"*^''^ l^f 1 '^^1 
unbound. [PB] partially bound and [IB] incorrectly bound. Thus [CB] Dji j) + [UB] + IPB] Opj) + [IB] <l.,ij) -TOpj, 

within 0.0001%. ........ t„, 

[0105] The model analytic for the estimated fluorescence activity V in a hixel h^^^ is given by 



45 t to (hybridization start time). 

[0106] The fluorescence stabilization section is given by 



^tH =H'k (/+6),5^0; 



[0107] And the rate of fluorescence stabilization is expressed via 



[0108] 



While fhfijj = ^ hoj) + 5 ) .S ^ 0: the system is assumed to be memoryless. i.e.. e is ergodic. 



[0109] The resulting dot spectrogram generated at step 206 is given by; 
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<I>(i.j) 

Preconditioning of The Dot Spectrogram 

' 101101 With reference to FIG. 5. the dot spectrogram *(i.i) is preconditioned by performing the following steps. First, 
at step 502. the dot spectrogram is refocused to yield a refocused dot spectrogram. Then, at step 504. a cross-^o^e^ 
lation convolution operation is perfom,ed to yield a correlated refocused dot spectrogram^ A local max.ma fitter 6 .s 
then applied at step 506 to the correlated refocused dot spectrogram to yield a maximized dot spectrogram. The max- 

10 Imized dot spectrogram is re-scaled at step 508 to yield a unifonnly re-scaled dot spectrogram. The unrfomily re-scaled 
dot spectrogram is then further re-scaled at step 510 by amplifying local edge hix^ boundaries o ""^^^^^^ 
scaled dot spectrogram to yield a globally re-scaled dot spectrogram denoted by* (/. fi. Also, amplrtude wandenngs 
within the globally re-scaled dot spectrogram are estimated at step 512 for use downstream, 
romi The steps of FIG. 5 result in a dot spectrogram tessellation (DST) operator. Purpose of DST Operator is to 

75 detemiine idealized ensemble boundaries for forcing downstream resonant action. Typically, in signal P'W'e^.ng af^ 
pSns.high pass or band pass spatial filtering is implemented toenhanceSNRin dot spectrogram matnx^^^^^^^ 

methods ;pply a combination of Laplacian or other edge detection filters apply to enhance signal from arrays cells w^ 
a h gher hSization concentration from those of the adjacent cells. These SNR enhancement methods however work 
only wtth Srtive or zero-SNR. Since SNR in general is negative in our case (ultra-low target DMA concentrations). 
20 these methods in effect amplify noise or further blur the hixelboundanes. ^, . , • r ^ „ 

01121 With the exemplary method, a new technique is used to amplify highly local morphological vana ions in the 
overall dot spectrogram matrix. A filter ^ is applied to the nonnalized amplitude dot spectrogram matrix *. for all com- 
binations +i+j+k where k is typically ranges from 0 to 2. 

25 Refocusing The Dot Spectrogram 

[0113] Refocusing of the dot spectrogram at step 502 to yield a refocused dot spectrogram is performed by deter- 
mining a locally averaged amplitude sutvarray represented by fl. then for each value of i. j and k. "^ -^"SJ^^ J^^ 
0 to 2: determining a local standard deviation o over a (2k*1 )X(2k*1 ) hixel neighborhood centered at (..j) and applying 
30 a filter ^ to the dot spectrogram *, for all combinations ±i±]±k. 
[0114] a is calculated by 



40 



45 



[0115] The calculation of n is independent of the oxel layout. From a precondittoning effectiveness standpoint, a 
best case design con^esponds to a completely random oxel layout in temis of 6-value Ibr adjacerrt oxels. The worst 
case corresponds to 6-value separation of 1 among adjacent oxels. The oxel. o(..j), . centered at (..j) compns^ of 
complementary oligonucleotides, corresponding to a mutation of interest over the set Z. l.oa. r. for averaging amplitude 
ranges from ±5 oxels to ±k oxels depending on (d-value mod 4^) separation, fl^z then is computed using 

50 When the averaging loci for multiple mutations of interest overlaps, then composite loci, f is bounded by a rectangle 
whose too left hand and bottom right hand comer coordinates are given by Bmin-k-1 .imax+k* 1] and linnax+K*^ .Imin ^-^ J 
where i^^. i„««. j^n and w correspond to the ordinate and abscissa Ibr the oxel mutattons with overlapping loci. 



'min' 'max* 

[0116] C given by: 



55 0(i j)'^ = 0(i j) . K , [ n(<I>(i j).k) /o(<I>(i,j),k)l V i.i 

10117] The preceding expression captures local ensemble deviation from dot spectrogram average. The design pa- 
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rameter k is computed offline based on specific bic^molecular signatures of interest The suffix / denotes iterative index 
for a cross-con-elation convolution operation to be applied after global refocusing. 

[01181 A uniform dot spectrogram reseating can be achieved by applying a. d(l,j) where a can be either a constant 
or a functional. This operator selectively enhances those hixels and ensemble boundaries whose intensity exceeds 
the local average by more than Kprespecified standard deviations. v^^^^^eon 
[01 19] K is a predetermined constant computed offline based on microarray fluorescence or chemiluminiscence sen- 
sitivlty. 



. Y loulloUgomoclcoiidcstaul \ 



[0120] The suffix / denotes an iterative index for a cross-conrelation convolution operation to be applied after global 
refocusing. 

Cross-con-elation of the Dot Spectrogram 

101211 The cross-correlation convolution operation of step 504 to yield a correlated refocused dot spectrogram is 
nerfonned by convolving the refocused dot spectrogram with an apriori-chosen restricted field obtained using a strongly 
dissipative dynamical system. Dissipative dynamical systems are those which define a fonA^ard regularizing flow in an 
adequate phase space containing an absorbing set. An absorbing set is a bounded set that attracts all bounded solu- 
tions in finite time at an exponential rate. Since we exploit a strongly dissipative system, the absorbing set is required 
to be unique compact set that is both positively and negatively invariant under the flow, it attracts all flows 
[0122] Conceptually, the post-hybridization hixel denotes a projection of the flow field that is absorbed by the coupling 
system Since we couple with an absorbing set. this stage yields significant SNR enhancement. 
[01231 In general, the convolving field can be constmcted using Kuramoto-Sivashinsky equation. 2D Navier-Stokes 
eouation or some forms of Reaction-Diffusion equations. The dot spectrogram suban^y around the oxel detecting 
mutation of interest can be cross-correlated with any dissipative dynamical system. In summary this step exploi s 
proven classes of mathematical system or exponential attractors. By definition exponential attractor is an exponentially 
contracting compact set with finite fractal dimension that is invariant under the forward flow. 

[0124] In the exemplary implementation, the following is performed to achieve the cross-conrelation. A synthetic 
system of the form is assumed, whose evolution equation is given by 

u, + Au + R(u) = 0 
u(0) = Uq 

where A is a positive self-adjoint operator with compact inverse and nonlinear term R satisfies the order match 
condition in step 306 in Figure 3. . 
[0125] The functional formulation of2D Navier-Stokes equations for incompressible fluid flow is chosen as the cou- 
pling system 
[0126] It is given by 

Ut + vAu + B{u, u) = f, 
u{0) = Uq 

where A = -PhA is the Stokes operator. B(u.u) stands for the non-linear terni (u. Au) projected to the underlying 
Hilbert space H. f is the volume force projected to the same Hilbert space and v is the viscosity temn. To meet the 
inoompressibility condition the following Hilbert space is assumed to set the initial value problem 
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H = I « : div u = 0 , Ju(x)dx = 0 , Uj},,,^ = u.| ,,„o. 4 ^ Uij 

and 

where Q Is the square [O.Ll X [O.L]. The mathennatical domain of the Stokes operator is given by 

D(A) = (Q)^ nV 

[01 27] Galerkin Approximation can be used to approximate the exponential attractor for the system 
[01281 The above system with periodic boundary conditions admits an exponential fractal attractor M in B whose 
dimensions can be estimated using 

dF(M) < c G^(log +1) 

where G is the Grashoff number given by and c is a constant that depends on the shape factor. denotes 
25 the positive eigenvalue of A. 

[0129] The rate of convergence of this system can be computed as well. 

[01 30] The exponential attractor is then coupled with the post-hybridization dot spectrogram subanray 
[0131] The exponential attractor is discretized over the grid that corresponds to the refocused amplitude subarray 
associated with a mutation and estimated above. 
30 [0132] The actual convolution of the two systems is then given by 
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[01 33] This process is computed for each oxel associated with the loci of a specific mutation of interest. 
40 Maximizing the Dot Spectrogram 

[0134] Step 506 for applying a local maxima filter ^ at to the correlated refocused dot spectrogram to yield a maxi- 
mized dot spectrogram is then implemented as follows. The local maxima filter is defined by 

n'=<S>'(U) x6(/j) 

wherein O(l.j) = 1 if OKij)' is a local maxima (or defines a local maxima ensemble with several neighboring hixels 
with same amplitude intensity); and wherein d<l j) = 0 othenwise. 



Re-scaling the Dot Spectrogram 



[0135] The maximized dot spectrogram is reseated at step 508 to yield a unifomily re-scaled dot spectrogram by 
applying an operator a . O(l.i) where a is either a predetemiined constant or a predetennined functional. 
55 [01 36] a is set so as to unifomily scale the convolved dot spectrogram to 1 /(I +exp(-MIN,NT/2)) to 1/(1 +exp(MAX,^,T' 
2)). Where [MIN,nt.MAX,nt) refers to the dynamic range of the microarray . . ^ u- i 

[01371 The unifomily re-scaled dot spectrogram is itself then rescaled at step 410 by amplifying local edge hixel 
boundaries of the unifomily re-scaled dot spectrogram to yield a globally re-scaled dot spectrogram. This is achieved 
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by 1) determining the zero mean amplitude for the uniformly re-scaled dot spectrogram; 2) applying a logarithmic 
rescaling function p around the zero mean amplitude; and 3) merging the local maxima into a single local maximum 
half way in between. ^ 
[0138] The logarithmic reseating function is generated by generating an expansion sequence of nonnegative num- 
bers and by generating an expanded dot spectrogram tessellation for O. 
[0139] The expansion sequence is generated as follows: 

iP)^ PO' P2^ - 
which is strictly decreasing until it reaches zero and thereafter is equal to zero by using: 



X'^ forn^O 
for/i<0 



20 



30 



[01401 The expanded dot spectrogram tessellation for O (which is represented by DST(O)) Is generated using: 

DST(<D) = p ( D{<l)i. 

wherein D{^,. ^) is a shortest possible discretized fluorescence amplitude separation between a pair of hixels 
25 (hi, Ak wherein <^y^ is a local maximum. u- i * 

[01411 The local maxima are merged into a single local maximum half way in between for downstream hixel-to- 
ensemble and ensemble-to-ensemble operations on hixel clusters using the sequence 



max{ « {(c) - K (m), 0 }; m = 0. 1 . 2,. 



where c is some positive integer constant. . ^ a 

[01421 Additional details are as follows, for any pair of dot spectrogram hixels, (D^ and <Dj, let the d.scretized fluores- 
cence amplitude separation between them be denoted by D(<t)i, (|)|). So, D(<t)i, <j)|) = 0 if and only if 

35 jdh. <b:| < C where C = Max(p mi))) - ^HP (^(i j)))/N and D((Di. <>j)=1 if and only if and <t»j are adjacent to each other 
on the normalized scale, but not identical, i.e.. |<|>i. ^ 2C. Let <Dk be a local maxima image for some other dot spec- 
trogram realization <S>(i,jT, ■ Define m,. <t>k) to be the shortest distance possible from hixel to some hixel which is 
a local maximum. Cleariy. D(4), 0^) = 0 ifand only if <Dk is a local maximum and D(<l»p ^) = I if and only if <Dk is not a local 
maximum but is adjacent to some local maximum. 

40 [0143] The an expansion sequence is thus defined to be any sequence of nonnegative numbers {p} = Po- Pa* 
(as mentioned above) which is strictly decreasing until it reaches zero and thereafter is equal to zero. For 
example, {p} = 15, 11, 5, 2. 1. 0, 0. 0, 0.... If O is a local maxima DST and {p} is an expansion sequence, then the 
expanded DST(O) is defined as 

^5 DST((0) = p ( D(t|>j, 0k)) 

[0144] K • Z ->Z on integers is useful for defining expansion sequences. A key property of this sequence is that given 
only the value of the sum s=p+q. where p and q are any numbers from this sequence, then we can always detemiine 
50 the values of the two original summands p and q. For example, choose the expansion sequence to be: 
31 15 7. 3. 1. 0, 0. 0. 0,... (generated by K(5-i): for i = 0.1.2, ...) 

and let DST. and DSTg be two DST realizations "expanded" using this sequence. This unique summation prop- 
erty ensures that a threshold value can be chosen for the DST summation DST^+DSTb so that the resulting realization 
retains only those local maxima in the two images which are within bounded intensity of n of each other. This expansion 
55 sequence preserves the individual local maxima from each of the two realizations, 

[0145] For downstream hixel-to-ensemble and ensemble-tOHonsemble operations on hixel clusters the local maxima 
are merged into a single local maximum halfway in between using the sequence using 
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max{ K{(c) - K(m), 0 ); m = 0. 1, 2,... 



where c Is some positive integer constant. 
[0146] In the alternative, a Dirichlet tessellation operator or a Delaunay triangulation operator are then appiiea lo 
perfomn gradient refocusing rather than steps 502 - 510. 

Estimation of Amplitude Wanderings 

[01471 Amplitude wanderings are estimated at step 412 within the enhanced dot spectrogram. The estimate is per- 
formed by applying a Palm Distribution operator to the globally re-scaled dot spectrogram to capture amplitude wan- 
derings and transitions at element, neighboring pair and local ensemble levels. The application of the Palm Distribution 
operator generates bounds that are used to accommodate degradation of hybridization overtime. Essentially the esti- 
mate exploits the use of generator, functions to capture stochastic variability in hybridization binding efficacy and draws 
upon results in stochastic integral geometry and geometric probability theory. 

101481 Geometric measures are constmcted to estimate and bound the amplitude wanderings to facilitate detection. 
In particular we seek a measure for each mutation-recognizer centered (MRC-) hixel that is invariant to local degra- 
dation. Measure which can be expressed by multiple integrals of the form 



where Z denotes the set of mutations of interest. In other words, we detemriine the function f(z) under the condition 
that m(z) should be invariant with respect to all dispersions ^. Also, up to a constant factor, this measure is the only 
one which is invariant under a group of motions in a plane. In principle, we derive deterministic analytical transfomiations 
on each MRC-hixel., that map en-or-elliptic dispersion bound defined on^^ (the two dimension Euclidean space - i.e.. 
oxel layout) onto measures defined onat. The dispersion bound is given by 

Log4 I' )■ 

[0149] Such a representation uniqueness facilitates the rapid decimation of the search space. It is implemented by 
instantiating a filter constmcted using measure-theoretic arguments. The transfomiation under consideration has its 
theoretical basis in the Palm Distribution Theory for point processes in Euclidean spaces, as well as in a new treatment 
in the prx)blem of probabnistic description of MRC-hixel dispersion generated by a geometrical processes. Latter is 
reduced to a calculation of intensities of point processes. Recall that a point process in some product space E X F is 
a collection of random realizations of that space represented as {{% fi). |e, G E. fj € F }. ^ ^ ,u 

r0150] The Palm distribution, n of a translation (T„) invariant, finite intensity point process in is defined to the 
conditional distribution of the process. Its importance is rooted in the fact that it provides a complete probabilistic 
description of a geometrical process. * 
[01 511 In the general forni, the Palm distribution can be expressed in temis of a Lebesgue factonzation of the form 

EpN*= A Ln X n 

[0152] Where n and A completely and uniquely detemnine the source distribution P of the translation invariant point 
process Also Bp N* denotes the first moment measure of the point process and U, is a probability measure, 
[01531 Thus we need to determine n and A which can uniquely encode the dispersion and amplitude wandering 
associated with the MRC-hixeL This is achieved by solving a set of equations involving Palm Distribution for each 
hybridization (i.e., mutation of interest). Each hybridization is treated as a manifestation of a stochastic point process 

111 

[01 541 In order to determine n and A we have implemented the following measure-theoretic filter 
Determination of A 

[01 551 using integral fomiulae constructed using the marginal density functions for the point spread associated with 
MRC-hixel(ij) 
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[01561 The oligonucleotide density per oxel Pmi'Jl PGR amplification protocol (c^), A^^^^^^"^*'^^ 

(T, J and imaging performance provide the continuous probability density function for amplrtude wandenng in the 

m-th MROhixel of interest. Let this distribution be given by (pm(ij). .Tim • ^"m)- 

[0157] The method requires a preset binding dispersion limit to be provided to compute A. The second moment to 
the function 

at SNR = 0 condition is used to provide the bound. 

Determination off n 

[0158] Obtained by solving the inverse problenn 

n = e*p 

where 

where and represent the normalized hybridization dispersion limits. These number are empirically plugged 
in. We choose 0.1 and 0.7 respectively to signify loss of 10% - 70% hybridization. 
[01 59] Also, e denotes the distribution of known point process. We use the form 

[01601^^^Um™t2y"the^fi^ (or enhanced) dot spectrogram generated by the step of FIG. 5 is repre- 

sented by: 

Resonant Interaction 

[0161] Referring again to Fig. 2, at step 210. the resonant interaction between the QEF and the preconditioned dot 
spectrogram is performed until a pre-setected degree of convergence is 

achieved. The resonance interaction proceeds by iteratively calculating a resonance equation until convergence is 

achieved. . , ♦ ^ 

[0162] More specifically, the following resonance equation is iteratively calculated. 

du 

[01631 In the above equation V(u) is the actual precondition, refocused MROhixel subarray. So a stable equilibrium 
Sie(miai)a.^ay)istJ^^ 

[o^jTn me forgoing, "a" is a constant « fluorescence decoherence timescale(s). State variable u corresponds to 
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nonstationary Markov random field (NS-MRF). Y(t) = Asin (wt ^ t) where t Is a small random P^^^^^^lf °^'3!!.nH^^^^ 
function (control parameter that influences convergence rates), "t" denotes integration '"^^f ^'^^^^^^^ 
the log (PSD maxima) ofthe ground tmth order function summed over all regions of mterest. More precisely. QN(t) s 
reveriberation projection (at some instant t prior to resonant convergence) for the coupled Ofs d.scussed alx)ve .n 

[OitS" *T?e^pflm generators are used as additive correction temis to the potential gradient to compensate for un- 
certainty and post hybridization decay. The actual dynamics is given by V(u, n, X). 
[01 66] The resonant iteration is terminated when 
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or when iteration counter t exceeds 1 03 (digital approximation to analog dynamics). 
101671 In this method the output readout from the resonance interaction is trivialized to deliver point-of-care (POC) 

20 diagnostics systems. Atter exponential SNR enhancement (i.e.. resonance) only those k MRC-hixels will have an am- 
plitude which are actually present in the unknown target sample. 

F0I68] In the event no resonance is achieved during step 210. a new QEF (coaesponding to a new mutation of 
interest) is selected and the method reinitializes to the original DS as computed in step 508 in FIG 5. This technique, 
referred to herein as a "software programmable QEF flush control melhodotogy". of reinitializing the dot spectrogram 

25 and using a new QEF leads to a new cycle through a signal sorter for rapid single- / multi-point gene^muta ion sorting. 
It is amenable to RF. electronic as well as optoelectronic QEF loading to microelectronics computing backplane micro- 
array readout and analysis backplane itself optoelectronically or electronically bonded to btoelectionic substrate. Fur- 
themiore the method can be implemented on offline miniature custom VLSI / palmtop / desktop setup. The QEF is also 
implementable in electronic, optoelectronic, bionic and ionic COTS/custom device physics In ^^f'^l^^^^;^ ^ 

30 designed such that complementary oligonucleotides for several mutations are spatially well spaced out «hat mul ip e 
QEF can be introduced to stimulate the system at the same time, then a logical "OR".ng of the resonance output is 
used to reach the detection decision. x^^i,^tt^,, 
101691 Also note that a chip fabricated to implement the exemplary method can wori< in two ways: a) seek all mu- 
tations of interest simultaneously: or b) seek all mutations of interest serially An advantage of (a) over (b) .s compu- 

35 tational speedup. An advantage of (b) over (a) is serendipity, i.e.. in method (a) only those mutations 

amplified that are detectton candWates. Everything else is likely suppressed or decimated. In (b) m"»iP'e -^sonan 
outputs can be accepted. By accepting multiple peaks for the OF the method can actually accept 2"-* orde . 3^ order 
entrained states, where order implies hamming distance to the mutation of interest in terms of base pair labels and 
locations. This can be used to theoretically accept entire families of denved mutations. 

40 roi701 Preferably, the resonance interaction is performed digitally by applying a matnx representative of the reso- 
nance equation to a matrix representative of the resonance stimulus in combination with a matrix representative of the 

lOITlT TffJfinal result of step 210 is a set of hixel locations wrtierein resonance has occurred Wentified by row and 
column number (i.j). 

Resonance Output Interpretation to Identify Diseases 

101721 Once resonances peak are observed in specific hixel locations of the preconditioned dot spectrogram at step 
210 the hixel addresses (k.i) of those locations are mapped at step 212 into the oligonucleotide table mentioned above 
(which contains the oligonucleotide sequences associated with hixel locations) to thereby identify the mutations, if any 
present in the sample being analyzed. This is a simple table look-up i«sultk,g In a direct readout of the mutations Fw 
a custom POC diagnostic sensor only those hixels which relate to mutaltons or expressed products of interest are 

mi731 At step 214. the mutations identified at step 212 are then mapped into the mutations table mentioned above 
which contains the diseases associated with the mutations) to thereby Wentify the diseases, if any, present in the 
sample being analyzed. This is also a simple table look-up resulting in a direct readout. Again, for a custom POC 

diagnostic sensor only those mutations which relate to diseases of interest are stored in the table. 

(0174) More specifically, at steps 21 2 and 214. the resonant output interaction is interpreted to yield a set of confimied 
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mutations as follows: ^.u -j- ^»..^i^*:ha 

MRC-hixel 2 = o(h(kJ)) where o is some hashing function or table look-up. The step will readout the digonucleotiae 
sequence from a table that has encoded the microan^y. Note that this mapping also operates to map the confirmed 
mutations to known diseases associated with the pre-selected set of mutations of interest to identify diseases, if any, 
indicated by the DNA sample (step 214). Note also that no probabilistic inferencing. exploitation of learning or nonlinear 
mapping is required to interpret the resonance output. Rather interpretation is rendered very straightforward thereby 
requiring only a low-cost hardware implementation with simple software to implement steps 212 and 214. So an im- 
portant aspect of the exemplary method is to accurately and robustly detect specific oligonudeottde sequences in the 
target sample. Subsequent association to understood genomics pathways is trivialized. 

Diagnostic Confirmation 

10175] Finally the diagnosis generated by step 214 is confirmed at step 216 by taking the identified diseases and 
solving in reverse for the associated QEF and then comparing the associated QEF's with ones expected for the mu- 
tations associated with the identified diseases to verify con^espondence and. if correspondence is not found, then a 
new set of mutations of interest are selected and all steps repeated. \ -tk- • 

[0176] In other words, this step maps detected mutations and expressed genes to a diagnostic assessment. This is 
a probabilistic or detemiinistic step, depending upon the genomics of the specific disease. It is represented as 

Diagnosis p fa >o - A/mutationj, mutationj mutation,) 
[0177] Where ^ is the specific diagnostic model associated with a medical condition of interest. 
[0178] Step 216 is a double check mechanism adopted to confimn multi-factorial diseases where the biochip encodes 
complex genomics. 

Aiternative Embodiments 

[0179] What has been described thus far is an implementation whereby a phase-space resonance interaction is 
performed. This implementation, which may also be referred to as a "mixed-mode phase shifted mode" is particularty 
effective for automatically extracting an entire dass of mutations that may be manifested in a hybndized element, in 
general the mixed-mode provides polymorphism in induced couplings for QEF design which delivers repeatability on 
analysis whereby mutation signatures of interest are simultaneously coupled to many base "dynamical systems vvith 
a single phase-embedding operator. Other resonance coupling interactions may be exploited as well. Other examples 
of couplings, other than phase-based, are "additive coupling mode" which provides further SNR enhancement and a 
"shunted input multiplicative coupling mode" which amplifies noise-to-noise couplings and leads to denvation ofbetter 
readout threshold for diagnostics dedsion making. Also, a combination of different resonance interactions can be 

Soi^betails regarding a related implementation may be found in U.S. Patent 6,142.681 entitled "Method and Ap- 
paratus for Interpreting Hybridized-Bioelectronic DNA Microarray Patterns Using Self Scaling Convergent Reverberan 
Dynamics". Details regarding an implementation directed to measuring viral loads may be found in U.S. Parent 
6,245.51 1 entitled "Method and Apparatus for Exponentially Convergent Therapy Effectiveness Monitonng Using DNA 
Microanray Based Viral Load Measurements". ^ ^ ... * »• 

[01 811 The exemplary embodiments have been primarily described with reference to flow charts illustrating pertinent 
features ofthe embodiments. Each method step may also represent a hardware or software component for performing 
the con-esponding step. It should be appredated that not all components or method steps of a complete implementation 
of a practical system are necessarily illustrated or described in detail. Rather, only those components or method steps 
necessary for a thorough understanding of the invention have been illustrated and described in detail. Actual imple- 
mentations may utilize more steps or components or fewer steps or components. 



Claims 

1. A method for analyzing an output pattern of an arrayed microdevice to identify mutations, if any. present in a 
biological sample applied to the arrayed microdevice said method comprising the steps of: 

generating a resonance pattem representative of resonances between a stimulus pattem assodated with a 
set of known mutations and the output pattem; and 

interpreting the resonance pattem to yield a set of confinmed mutations by comparing resonances found therein 
with predetemriined resonances expected for the set of mutations. 
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2 The method of daim 1 wherein the output pattern is derived from an arrayed micfode>«ee ►'^^^"S,? P'"'^''^ °! 
each having a set of substantially identical immobilized oligonucleotides with the oligonucleotides of each cell 
being substantially unique from the oligonucleotides of every other cell. 

3. The method of daim 1 wherein the step of generating a resonance pattem indudes the step of iteratlvely processing 
the output pattem by performing a convergent reverberation to yield a resonance pattern ^^P^f ^"'^Jrf ° ^^^^^^ 
nances^between a predetemiined set of selected Quantum Expressor Functions and the output pattem until a 
predetermined degree of convergence is adiieved between me resonances found in the resonance pattem and 
resonances expected for the set of mutations. 

4. The method of daim 2 further induding the steps of: 

selecting a sub-set of the mutations for analysis; r.- _Ki-.4»= 

generating preconditioner transfomis based upon the selected sub-set of mutations by inverse Dindilet tes- 
sellation on the arrayed microdevice structure; and . ^ ,. n 
selecting a sub-set of nonlinear Quantum Expressor Functions from a set of predetemiined nonlinear Quantum 
Expressor Functions based upon the selected sub-set of mutations. 

5. The method of daim I further induding the step of modifying the output pattem to generate a modified output 
pattem by perfomiing one or more of the steps of: 

refocusing the output pattem to yield a re-focused output pattem; 

cross-conrelating the re-focussed output pattem; . . j . 

applying a local maxima filter to the correlated re-focussed output pattem to yield a maximized output pattem, 

re-scaling the maximized output pattem to yield a unifomily re-scaled output pattem; and 

re-scaling the uniformly re-scaled output pattem to amplifying local boundaries therein to yield a globally re- 

scaled output pattem. 

6. The method of daim 5 wherein the step of re-seating the maximized output pattem is performed by aPP^"!^" 
operator to selectively enhance those hotels and ensemble boundaries whose intensity exceeds the local average 
by more than a prespecified number of standard deviations. 

7 The method of daim 5 wherein the step of reseating the unifomily re-scaled output pattem is performed by applying 
a logarithmic reseating function around zero mean amplitudes within the output pattem to amplify local edge bound- 
aries of the output pattem. 

8. The method of daim 1 further including the step of 

estimating amplitude wanderings within the globally re-scaled output pattem. 

9. The method of daim 8 wherein the step of estimating amplitude wanderings is perfomned by applying spectrogram 
ensemble invariance property generators constructed using a Palm distribution. 

1 0. The method of daim I wherein the step of interpreting the resonance pattem to a yield a set of confirmed ^^tions 
by comparing resonances found therein wHh predetenmined resonances expected for the set of mutations indudes 
the step of 

mapping the confirmed mutations to knovwi diseases associated with the pre-selected set of known mutations 
to identify diseases, if any, indicated by the biological sample. 

11 The method of daim 1 further induding the step of perfbrming a diagnostic confinnation by taking the identified 
■ diseases and solving in reverse for the assodated Quantum Expressor Functions and then ~;"P^""9f:°^f 
turn Expressor Functions with ones expected for the mutations assodated with the identified disease to venfy 
correspondence and, if no correspondence is found, selecting a new sublet of known mutations and repeating 

all steps. 

12. The method of step 1 further induding the initial steps of generating the set of nonlinear Quantum Expressor 
Functions by: 



24 



EP 1 145 181 B1 

inputting a set of mutation signatures representative of the pre-selecled set of known mutations; 

inputting a representation of an arrayed microdevice oligonucleotide pattern layout for the arrayed microdevice 

from which the output pattem is generated; 

generating a set of resonant interaction parameters representative of mutation pattem interactions between 
elements of the arrayed microdevice including interactions from a group including element-to-element inter- 
actions, element-to-ensemble interactions, ensemble-to-element interactions, and ensemble-to-ensemble in- 
teractions; and . * « 
generating the set of nonlinear Quantum Expressor Functions from the set of resonant interaction patterns. 

13 The method of claim 12 wherein the step of generating the set of nonlinear Quantum Expressor Functions com- 
prises the steps of matching a power spectral density (PSD) amplitude of a coded mutation signature, correspond- 
ing to the pre-seleeted mutation set of interest, to that of Hamiltonian system so that stochastic and detemiinistic 
time scales match, and the time scales couple back to noise statistics and degree of asymmetry. 

14. The method of step 1 further including the initial steps of generating the set of nonlinear Quantum Expressor 
Functions by: 

calculating values representative of a pre-selected Hamiltonian function; 
calculating hamrionic amplitudes for the Hamiltonian function; 
generating an order function from the Hamiltonian; function, 
measuring entrainment states of the order function; and 

modulating the order function using the entrainment states to yield the Quantum Expressor Function. 

15. The method of step 14 wherein the Hamiltonian is a Spin-Boson Hamiltonian. 

16. The method of claim 14 wherein the step of generating the order function includes the step of approximating a 

Daido Integral. 

17. The method of claim 1 wherein the output pattem is a quantized output pattem. 

18. The method of claim 17 wherein the quantized output pattem is a dot spectrogram. 

19. The method of claim 1 wherein the biological sample selected from a group consisting ofa DNA. RNA. protein, 
peptide-nudeic acid (PNA) and targeted nucleic amplification (TNA) samples. 

20. A system for analyzing an output pattem of an arrayed microdevice to identify mutations, if any. present in a 
biological sample applied to the arrayed microdevice. said system comprising: 

a resonance pattem generator operative to generate a resonance pattem representative of resonances be- 
tween a stimulus pattem associated with a set of known mutations and the output pattem; and 
a resonance pattem interpretation unit operative to interpret the resonance pattem to yield a set of confirmed 
mutations by comparing resonances found therein with predetemnined resonances expected for the selected 
set of mutations. 

21 The system of claim 29 wherein the output pattem is derived from an an-ayed microdevice having a plurality of 
cells each having a set of substantially identical immobilized oligonucleotides with the oligonucleotides of each 
cell being substantially unique from the oligonucleotides of every other cell. 

22 The system of daim 20 wherein resonance pattem generator includes a resonant interaction unit operative to 
iteratively processing the output pattem by perfomiing a convergent reverberation to yield a resonance pattem 
representative of resonances between a predetemiined set of selected Quantum Expressor Functions and the 
output pattem until a predetermined degree of convergence is achieved between the resonances found in the 
resonance pattem and resonances expected for the set of mutations. 

23. The system of claim 20 further including a image preconditioning unit operative to perfomi one or more of the 
functions of: 

refocusing the output pattem to yield a re-focused output pattem; 
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cross-correlating the re-focussed output pattern; . . ^ » *«^»^r«. 

applying a local maxima filter to the correlated re-focussed output pattern to yield a maxfmized output pattern, 
re-scaling the maximized output pattern to yield a uniformly re-scaied output pattern; and , ^ „ _ 

re-scaling the unifonnly re-scaled output pattern to amplifying local boundaries therein to yield a globally re- 
scaled output pattern. 

24 The system of daim 20 wherein the resonance pattern interpretation unit is operative to map the confin^ed mu- 
tationrto known diseases associated with the preselected set of known mutations to Klentify diseases, if any. 
indicated by the btological sample. 

25. The system of step 20 further including a Quantum Expressor Function generator operative to: 

input a set of mutation signatures representative of the set of known mutations; . kv-k ».« 

input a representation of a arrayed mlcrodevice pattern layout for the arrayed mioodevice from which the 

renerTaTet'XTna'n^^ parameters representative of mutation pan^ inte^ctions t^tween 

elements of the arrayed microdevice Including interactions from a group including ^l^f^" ""^^l"! "* "1*^ 
acttons. element-to^nsemble interactions, ensembto-to-element interactions, and ensemble-to-ensemble in- 

gen2ateVe"set of nonlinear Quantum Expressor Functions from the set of resonant interaction patterns. 

26. The system of daim 20 wherein the output pattern is a quantized output pattern. 

27. The system of daim 26 wherein the quantized output pattern is a dot spectrogram. 

28. The system of daim 26 wherein the biological sample is selected from a group consisting of a DNA. RNA. protein, 
peptlde-nudeic acid (PNA) and targeted nudeic amplification (TNA) samples. 

29. The method according to one of the daims 1 to 1 9 wherein an output pattern of a biodiip is used as the arrayed 
mlcrodevice output pattern. 

30. The system according to one of the daims 20 to 28 for analyzing an output pattern of a btodnip as the arrayed 
microdevice output pattem. 



Verfahren zum Analysieren eines Ausgabemusters einer Mikroarrayvorridntung, urn in einer auf die Mikroarray^^^^^ 
richlung aufgebrachten biologisdien Probe, wenn Qberhaupt. vorhandene Mutationen zu identifizieren. wobe. das 



Patentanspriiche 
1. 

Verfahren die Schritte umfasst: 

Erzeugen eines Resonanzmusters. das Resonanzen zwischen einer Gruppe bekannter Mutationen zugeord- 
neten Stimulusmuster und dem Ausgabemuster angibt. und ^ . .fno- 

interpretieren des Resonanzmusters, um eine Gaippe bestatigter Mutationen zu ert^atten. ^^^^^J^^^^^ 
fundene Resonanzen mit fOr die Gmppe von Mutationen envarteten. vorbestimmten Resonanzen verghchen 
warden. 

2 Verfahren nadi Ansprudi 1. bei dem das Ausgabemuster von einer Mikroarrayvorriditung erhalten wird. die eine 
MerrLh?Zel^VulWeist.d ejeweils eine Gruppe^ 

i;!Ss^.>!«^MirOlionukleoidejederZellegegenQberdenOlionukleotideniederander^ 
eindeutig sind. 

3 Verfahren nach Anspruch 1. bei dem der Schritt. ein Resonanzmuster zu erzeugen, den Schritt umfasst, dass 
AusgaSmusL iterativ zu verart,eiten. indem ein konvergentes Nadihallen durchgefuhrt wird. um e,n Reson^- 
muster zu erhalten, das Resonanzen zwis*en einer vorbestimmten Gmppe '^^^^^^'^^^'''''^^Z^ 
Funktlonen und dem Ausgabemuster angibt. bis ein vorbestimmter Konvergenzgrad zv«sd,en den 
nanzmuster aufgefundenen Resonanzen und fur die Gmppe von Mutationen envarteten Resonanzen erreidrt w,«J. 
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4. Verfahren nach Anspmch 2, femer die Schritte umfassend: 

Atiswahlen einer Unterqmppe der Mutationen zur Analyse, , . , ^. ^ 

S^v^'^dingungstransformiertenbasierendaufd 

inupnsA Dirichlet-Tesseliemng auf der Mikroarrayvorrichtung, und 

SZ^r^^uJ-E^^-^^r^^^ basierend auf der ausgewShlten Untergruppe Mutal.onen. 

5. verfahren nach Anspmch 1 . femer den Schritt umfassend. das ^^s^.^f ^^"^"^f '^^^^ 
ziertes Ausgabemuster zu erzeugen. indem einer Oder mehrere der Schntte durchgefuhrt warden. 

Refokussieren des Ausgabemusters. urn ein refokussiertes Ausgabemuster zu erhallen. 

Se'^^r^TstLT^^^^^^ urn eIn g.eichm..>g reskaUer.es Ausgabemuster zu erha.ten. 

Retkalieren des gleichmaBig reskalierten Ausgabemusters. um lokale Grenzen darinnen zu verstaiken. urn 
ein umfassend reskaliertes Ausgabemuster zu erhalten. 

6 Verfahren nach Anspruch 5, bei dem der Schritt. das maximierte Ausgabemuster zu reskalieren. durchgefahrt 
^r? in'er^ ein Sator arigewendet wird. um diejenigen Hixel und Ensemblegrenzen selekt.v zu ve.3tarf.en 
Sn intensitat den lokalen Mittelwert um mehr als eine vorbestimmte Anzahl von Standardabweichungen uber- 
schreitet. 

7 verfahren nach Anspruch 5. bei dem der Schritt. das gleichmSBig reskalierte Ausgabemuster 2^ T^^'^'^u^' 
glSwirXLemeinel^ 

bemuster angewendet wird. um lokale Randgrenzen des Ausgabemusters zu versfari<en. 
8. Verfahren nach Anspruch 1 . femer den Schritt umfassend. 

Amplitudenwanderungen In dem umfassend reskalierten Ausgabemuster abzuschStzen. 

9 Verfahren nach Anspruch 8. bei dem der Schritt. Amplitudenwandeoingen abzuschatzen, durchgefuhrtwird. indem 
untefve[!^nl^^ Palm-Verteilung aufgebaute Erzeugende fur Invarianzcharakterist.ka des Spektrogram- 
mensembles angewendet werden, 

10 Verfahren nach Anspruch 1, bei dem der Schritt. das Resonanzmuster zu inlerpretieren. um eine Gmppe besta- 
ligter Mutationen zu ertnalten. indem darin aufgefundene Resonanzen mit fur d.e Gruppe von Mutationen erwar- 
teten, vortjestlmmten Resonanzen vergllchen werden, den Schntt umfasst. 

die bestatlgten Mutationen auf bekannte Krankheiten abzubilden. die der ^"^^ 

kannter Mutationen zugeordnet sind. um von der biologischen Probe, wenn ubertiaupt. angegebene Krank- 
heiten zu identifizieren. 

1 1 Verfahren nach Anspruch 1 , femer den Schritt umfassend, eIne diagnostische Bestatigung durchzufOhren. indem 
d^SSe^n^rkhelten herangezogen und die zugeordne^ 

f^Rtehtung geiast und dann diejenigen Quantum-Expressor-Funktionen mit denjen.gen verghchen werden^ d.e 
Z SStionen eo«artet we«len. die der Identifizlerten Krankheit zugeordnet sind, um e,ne Ubere,nst.mmung 
zuO^lSrund eine neueuntergruppebekannter Mutationen ausgewahit und alleS^ 

wrenn keine Obereinstimmung festgestellt vwrd. 

12. verfahren nach Anspmch 1 , femer den anfanglfchen Schritt umfassend. die Gruppe nichi linearer Quantum-Ex- 
pressor-Funktionen zu erzeugen durch: 

Eingeben einer Gmppe von Mutetionssignaturen. die die vorab ausgewShlte Gmppe bekannter Mutetionen 
EilS*^ einer Wiedergabe einer Mikroarrayvorrichtungsolionukleodidmusterauslegung fOr die Mikroar- 
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rayvorrichtung, aus der das Ausgabemuster erzeugt wird, 

Erzeugen einer Gajppe von Resonanzwechselwirkungsparametem, die Mutationsmusterwechselw>r1cungen 
zwischen Elementen der Mikroarrayvomchtung. einschlieRlich von Wechselwirkungen einer Gmppe angeberv 
die Elemenl-Element-Wechselwirkungen. Element-Ensemble-Wechselv^irKungen. Ensemble-Element-Wech- 
selwirkungen und Ensemble-Ensembie-Wechselwirkungen umfasst. und 

Erzeugen der Gruppe nicht linearer Quantum-Expressor-Funktionen aus der Gruppe von Resonanzwechsel- 
wirkungsmustem. 

13 Verfahren nach Anspruch 12. bei dem der Schritt. die Gruppe nicht linearer Quantum-Expressor-Funktionen zu 
' erzeugen. die Schritte umfasst. eine Leistungsspektra!dichte-(PSD)-Amplitude einer kodierten MutatK)nss.gnatur. 
die der im Interesse stehenden. vorab ausgewShlten Mutationsgruppe entspricht. auf die eines Hamilton-Systems 
abzubilden. so dass stochastische und deterministlsChe ZeitmaBstabe ubereinstimmen und die Zeitmalist&be mit 
Rauschstatistiken und einem Symmetriemaft rOckgekoppelt sind. 

14. Verfahren nach Anspmch 1 . femer die anfanglichen Schritte umfassend. die Gruppe nicht linearer Quantum-Ex- 
pressor-Funktionen zu erzeugen durch: 

Berechnen von Werten, die eine vorab ausgewahlte Hamilton-Funktion angeben, 
Berechnen von hanmonischen Amplituden fur die Hamilton-Funktion. 
Erzeugen einer Ordnungsfunktion aus der Hamilton-Funktion, 

Messen von Frequenzmitnahmezustanden der Ordnungsfunktion. und ^ * c 

Modulieren der Ordnungsfunktion unter Venwendung der Frequenzmitnahmezustande, um die Quantum-Ex- 
pressor-Funktion zu erhalten. 

15. Verfahren nach Anspmch 14. bei dem die Hamiltonianische eine Spin-Boson-Hamiltonianische ist. 

16. Verfahren nach Anspruch 14, bei dem der Schritt. die Ordnungsfunktion zu erzeugen. den Schritt umfasst. ein 
Daido-lntegral abzuschStzen. 

17. Verfahren nach Anspmch 1 . bei dem das Ausgabemuster ein quantisiertes Ausgabemuster ist. 

18. Verfahren nach Anspmch 17. bei dem das quantisierte Ausgabemuster ein Punklspektrogramm ist. 

19. Verfahren nach Anspmch 1 . bei dem die biologische Probe aus einer Gmppe ausgewahit ist. die aus DNA-. RNA-. 
Protein-, Peptid-Nukleinsaure-(PNA)- und Zielzenkemverstarkungs-(TNA)-Proben besteht. 

20 System zum Analysieren eines Ausgabemusters einer Mikroarrayvonichtung. um in einer auf die Mikroarrayvor- 
■ richtung aufgebrachten biologischen Probe, wenn uberhaupt. vorhandene Mutationen zu identifizieren. wobei das 
System umfasst: 

einen Resonanzmustergenerator. der betriebsfahig ist. ein Resonanzmuster zu erzeugen. das Resonanzen 
zwischen einer Gmppe bekannter Mutationen zugeordneten Stimulusmuster und dem Ausgabemuster angibt. 

eilfe Resonanzmusterinterpretationseinheit. der betriebsfahig Ist. das Resonanzmuster zu interpretieren. um 
eine Gmppe bestatigter Mutationen zu ertialten. indem darin aufgefundene Resonanzen mit fQr die Gmppe 
von Mutationen enwarteten, vorbestimmten Resonanzen verglichen werden. 

21. System nach Anspmch 20. bei dem das Ausgabemuster von einer Mikroarrayvorrichtung erhalten wird. die eine 
Mehrzahl Zellen aufweist. die jeweils eine Gmppe im Wesentlichen identischer. immobilisierter Olionukleotide auf- 
weisen, wobei die Olionukleotide jeder Zelle gegenuber den Olionukleotiden jeder anderen Zelle im Wesentlichen 
eindeutig sind. 

22 System nach Anspmch 20. bei dem der Resonanzmustergenerator eine Resonanzwechselwirkungseinheit um- 
fasst die betriebsfahig ist, das Ausgabemuster iterativ zu verarbeiten, tndem ein konvergentes Nachhallen durch- 
gefuhrt wird, um ein Resonanzmuster zu erreichen. das Resonanzen zwischen einer vorbestimmten Gmppe aus- 
gewahlter Quantum-Expressor-Funktionen und dem Ausgabemuster angibt. bis ein vorbestimmter Konverpenz- 
grad zwischen den Resonanzen, die in dem Resonanzmuster aufgefunden werden, und Resonanzen erreicht wird, 
die fur die Gmppe Mutationen erwartet werden. 
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23. System nach Anspmch 20. femer eine Bildverarbeitungseinheit umfassend. die betriebsfShig ist. eine oder meh- 
rere der Funktionen auszufuhren: 

Refokussieren des Ausgabemusters. urn ein refokussiertes Ausgabemuster zu ertialten. 
Kreuzkorrelieren des refokussierten Ausgabemusters, , 
Anwenden eines Filters fur lokale ISIaxima auf das korrelierte refokussierte Ausgabemuster. urn ein maximier- 
tes Ausgabemuster zuertialten, 

Reskalieren des maximierten Ausgabemusters. urn ein gleichmafiig reskaliertes Ausgabemuster zu erhalten, 

Reskalieren des gleichmaBig reskalierten Ausgabemusters. urn lokale Grenzen darinnen zu verstarken. urn 
ein umfassend reskaliertes Ausgabemuster zu ertialten. 

24 System nach Anspmch 20, bei dem die Resonanzmusterinterpretatlonseinheit betriebsfShig ist. die bestatigten 
* Mutationen auf bekannte Krankhelten abzubilden, die der vorab ausgewahlten Gruppe bekannter Mutationen zu- 
geordnet sind. um von der biologischen Probe, wenn uberhaupt. angegebene Krankheiten zu identifizieren. 

25. System nach Anspruch 20, ferner einen Quantum-Expressor-Funktionsgenerator umfassend. der betriebsfahig 
ist, um: 

eine Gmppe von Mutationssignaturen einzugeben, die die Gruppe bekannter Mutationen angeben. 

eine Wiedergabe einer Mikroan-ayvorrichtungsolionukleodidmusterauslegung fOr die Mikroanayvomchtung 

einzugeben, aus der das Ausgabemuster erzeugt wird, 

eine Gmppe von Resonanzwechselwirkungsparametem zu erzeugen. die Mutationsmustenwechselwirkungen 
zwischen Elementen der Mikrcarrayvorrichtung. einschlieBlich von Wechselwirkungen einer Gmppe angeben. 
die Element-Element-Wechselwirkungen. Element-Ensemble-Wechselwirkungen. Ensemble-Element-Wech- 
selwirkungen und Ensemble-Ensemble-Wechselwirkungen umfasst, und 

die Gmppe nicht linearer Quantum-Expressor-Funktionen aus der Gmppe von Resonanzwechselwirkungs- 
mustern zu erzeugen. 

26. System nach Anspmch 20, bei dem das Ausgabemuster ein quantisiertes Ausgabemuster ist. 

27. System nach Anspmch 26. bei dem das quantisierte Ausgabemuster ein Punktspektrogramm ist. 

28. Verfahren nach Anspmch 26. bei dem die biologische Probe aus einer Gmppe ausgewahit ist, die aus DNA-. RNA-. 
Protein-, Peptid-Nukleinsaure-(PNA)- und Zielzellkemverstarkungs-(TNA)-Proben besteht. 

29. Verfahren nach einem der Anspriiche 1 bis 1 9. bei dem das Ausgabemuster eines Biochips als das Ausgabemuster 
einer Mikrcarrayvorrichtung venwendet wird. 

30. System nach einem der Anspruche 20 bis 28 zum Analysieren eines Ausgabemusters eines Biochips als das 
Ausgabemuster einer Mikroarrayvorrichtung. 



Revendications 



1 Proced6 d'analyse d'un motif de sortie d'un microdispositif ordonn6 pour identifier des mutations. !e cas 6ch6ant. 
pr6sentes dans un fechantillon biologique appliqu6 au microdispositif ordonn6. ledit precede comprenant les 6tapes 
consistent ^ : 

g6n6rerun motif de resonance repr6sentatif de r6sonances entre un motif de stimulus associ6 d un ensemble 
de mutations connues et le motif de sortie ; » • a 

interpreter le motif de resonance pour prt>duire un ensemble de mutations confirm6es en comparant les re- 
sonances d6couvertes dans celui-ci a des resonances pr6d6temiln6es attendues pour rensemble de muta- 
tions. 

2 Proc6d6 de la revendication 1 , dans lequel le motif de sortie est issu d'un microdispositif ordonn6 ayant une plurality 
de cellules ayant chacune un ensemble d'oligonucleotides quasiment identiques immobilises avec les oligonu- 
cl6otides de cheque cellule 6tant quasiment uniques par rapport aux oligonucleotides de chaque autre cellule. 
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3. PfOc6d6 de la revendicalion 1 . dans lequel r6tapc de production d'un motif de ^o^"** «^J^.2« 
sistent d trailer oar it6ration le motif de sortie en ex6ajtant une reveibSration convergente afin de produire un motif 
LS^nln^ S^^^^^^^^^ de resonances entre un ensembie pr.d.tem.n. de foncttons ^ Jl'P--^^^^^^^ 
s61ecUonn6es et le motif de sortie jusqu'd ce qu'un degre pr6d6lemi,ne de convergence ^"^^^^^' ^^^^ 
scnances d6couvertes dans le motif de resonance et des resonances attendues pour I ensemble de mutations. 

4. Proc6d6 de la revendication 2 comportant en outre les 6tapes consistant d : 

saectionnerunsous-ensemWede mutations pour analyse ; . . „„o .Aiowinnn* nar 

generer des transfbm,6es de pr6condHionnement ii partir du sous-ensemble de mutatwns s6lectionn6 par 
pavage de Dirichlet inverse sur la structure du microdisposrtif ordonne : et „„,„„hio h« 

selectionner un sous-ensemWe de fonctions d'expression quantique non lin^aires ^ Pa^-r d un ensemb^de 
fonctions d'expression quantique non Iin6aires pr6d6temiin6es sur la base du sous-ensemble de mutations 
s6lectionne. 

5 Precede de la revendication 1 comportant en outre Tetape consistant ^ modmer le motif de sortie afin de g6n6rer 
un motif de sortie modifi6 en executant une ou plusieurs des 6tapes consistant d : 

refocaliser le motif de sortie pour produire un motif de sortie refocalis6 ; 

20 effectuer une conr6lation croisee du motif de sortie refocalise ; 

appliquer un filtre de maxima locaux au motif de sortie refocalis6 correle pour produ.re un mot.f de sortie 

Illfeduire^e motif de sortie maximise ,. . 

rtduire le motif de sortie unifom^ement r6duit ^ des limites locales d'amplification a rint^rieur de celu-c. pour 
25 produire un motif de sortie globalement rSduit. 

6 Proc6de de la revendication 5. dans lequel I'etape de reduction du motif de sortie maximis6 est executee en 
a'Sant un oS'ateur pour majorer seleCivement les hixels et limites d'ensemble dont Hntensrte depasse la 
moyenne locale de plus d'un nombre pre-specifi6 d'6carts types. 

7. Precede de la revendication 5, dans lequel I'etape de reduction du motif de sortie ^'^'^'^^'l^J'''^*^^^ 
en appliquant une fonction de reduction logarithmique autour d'amplitudes de moyenne nulle d I .nt6neur du motif 
de sortie afin rfamplifier des limites de pourtour local du motif de sortie. 

8. Procede de la revendication 1 comportant en outre I'etape consistant d estimer des variations d'amplitude 6 I'in- 
terieur du motif de sortie globalement reduit. 

9 Precede de la revendication 8. dans lequel retape d'estimation de variations d'amplitude est executee en appl'^ 
qlTant des generateurs de proprietes d'invariance de I'ensemble des speclrogrammes constmits a I a«le d une 

distribution de Palm. 

10. Precede de la revendication 1. dans lequel I'etape d'interpretation du motif de '^'^'^^'J^^^^^^^ 
ble de mutations confim,6es en comparant des resonances decouvertes dans celu.-a d des ^O"^"^ P^J" 
terminees attendues pour I'ensemble de mutations comporte I'etape consistant d cartographier les '""tafcons con- 
Ss ^r reconnoitre des maladies associees ^ I'ensemble preseiectionne de mutations connues afin d iden- 
tifier des maladies, le cas ec»i6ant. indiquees par l'6chantillon biologique. 

11 Precede de la revendication 1 comportant en outre I'etape consistant ^i executer une confirmation de diagnosfc 
^S^siderantTes maladies identifiees et en resolvant a I'envers les fonctions d'expression ^"-^^'^^-^^^^ 
puis en comparant ces fonctions d'expression quantique ^ celles attendues pour '^^^^'^^^^^ 
maladie identmee afin de verifier une correspondance et. si aucune correspondence n est d6couverte. s6lect«nner 
un nouveau sous-ensemble de mutations connues et r6p6ter toutes les etapes. 

12. Precede de la revendication 1 comportant en outre I'etape initiate consistant a g6n6rer I'ensemble de fonctions 
55 d'expression quantique non lineaires en : 

introduisant un ensemble de signatures de mutations representatives de Pensembte pres6lectionn6 de muta- 
tions connues ; 
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teraction r^sonante. 

Bque non Iin6aires comprend I'^tape consistant S faire ""''„'''"^Se de Lte^^^ pr6s6lectionn6 

14. P«K*d6 de la revendication 1 comportant en outre .'..ape initia.e consistant . g.n^rer rensemble de fbocBons 
d'expression quantique non linSaires en : 

calculant des valeurs representatives d'une fonction hamiltonlenne pr6s6lectionn6e ; 
calculant des amplitudes hanmoniques pour la fonction tiamiltonienne ; 
g6n6rant une fonction d'oidre d partir de la fonction hamiltonienne ; 

TdS^r— ti^^^^^^^^^^^ 

15. Proc6d6 de la revendication 14. dans lequel I'hamiltonien est un hamiltonien de spin-boson. 

16. Pmc^.delarevendication14.danslequeir6tapedep«,duationdelafonctiond-ordrecomporter^^^^ 
d approximer une integrate de Daido. 

17. Proc6d6 de la revendication 1 . dans lequel le motif de sortte est un motif de sortie quantifi6. 

18 P«c6d6 de la revendication 17. dans lequel le motif de sortte quantifi6 est un spectrogramme de points. 

ungen^rateurde motif de resonance operantpourgenererun motif deresonan.^^^^ 
pr6d6termin6es attendues pour I'ensemble de mutations selectionnS. 

cellule. 

de mutations. 

23. Systemedeterevendication 20 comportant en outre une unite depreconditionnementtfimage operant pour 
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cuter une ou plusieurs des fonctions suivantes : 

refocaliser le motif de sortie pour produire un motif de sortie refocalis6 ; 

produire un motif de sortie globalement r6duit. 

25. Systftme de la «vendication 20 c»mportant en outre un g6n6rateur de foncUon d'expression quantique operant 
pour: 

':S::^°Z::Z.T.^^:Te^L^^ .u^m^e non nn^ires . partir de .-ensemble de mCHs d.nte- 
raction resonante. 

26. Systeme de la revendication 20. dans lequel le motif de sortie est un motif de sortie quantifife. 

27 SystSme de la revendication 26. dans lequel le motif de sortie quan«fi6 est un spectrogramme de points. 

2,.p™cM6»lo.^ruo.<»»v»dle.«oo.1S19.d=n.l«,uiu"mo«ia.»«aWpoc.4ADN.aut.eS^^ 

le motif de sortie du microdispositif ordonn6. 
30. systems se.on Tune des revendioations 20 a 28 pour analyser un motif de sortie d'une puce ^ ADN comme le 
motif de sortie du microdispositif ordonn6. 
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FIGURE 2 
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